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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1366 

HEAT TRANSFER BY FREE CONVECTION FROM HORIZONTAL 
CYLINDERS IN DIATOMIC GASES* 

By R . Hermann 


I. DETERMINATION OF HEAT -TRANSFER LAW FOR HORIZONTAL CYLINDER FROM 
TESTS WITH PARTICULAR ACCOUNT TAKEN OF THE 
TEMPERATURE CHARACTERISTIC Te^ 

1. Introductory Remarks 

The case of the horizontal cylinder is of particular importance in 
the study of heat transfer hy free convection for the following reasons: 
in the first place^ next to the rectangular plate it represents the 
simplest two-dimensional case; and second^ a very wide range of measure- 
ments is possible^ from the finest electrically heated glow lamp wires 
to pipes heated hy liquids or gases flowing through them. 

To investigate free-convection flow from the point of view of 
similarity considerations, it is convenient to consider the case of 
small temperature differences between the heated body and the surround- 
ings; in this case all the properties of the medium, even the density 
(ref.’l, p. 429), in the entire temperature field may be assumed 
constant. The case of large temperature differences, for which the 
variation in properties over the temperature range can no longer be 


’’Warmeiibergang bei freier Stromung am wagrechten Zylinder in 
zweiatomigen Gasen." VDI Forschungsheft, No. 379, 1936, pp. 1-24. 

^Accepted as dissertation by the Technical High School at Aachen 
with the approval of Prof. Dr.-Ing. C. Wieselsberger and Prof. Dr. 

W. Muller. The tests were carried out in the division for applied 
mechanics and thermodynamics of the Physical Institute at the University 
of Leipzig. The author takes this occasion to thank Prof. Dr. L. 
Schiller for the interest which he took in this work and for his valua- 
ble advice. 
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neglected, is treated separately. As first pointed out by W. Nusselt 
(ref. 2), it follows from the differential equations of free convection, 
together with the uniquely determining boundary values (ref. 1, p. 428) 
for the case of small temperature differences, that the nondimensional 
heat-transfer parameter Nu (Nusselt nvimber^) is determined by the lift 
coefficient Gr (Grashof number) and the characteristic Pr (Prandtl 
n\imber) for the molecular constitution of the gas^: 

Nu = F(Gr, Pr) (l) 


where 


Nu = od/X, Gr = d^ g P 0/v2 Pr = v/a (2) 

a is the heat-transfer coefficient (cal/cm^) (sec) (°C) , d is the cylin- 
der diameter (cm), g is acceleration of gravity (cm/sec^), X is the 
heat conductivity (cal/(cm) (sec) (°C) ) , v is the kinematic viscosity 
(cm^/sec), a is the temperature conductivity (cm^/sec), P is the coef- 
ficient of expansion (°C“^), and 0 is the temperature difference (°C) 
between body (t^) and medixzm (t„) . 

The case of large temperature difference between the body and the 
siirroundings, for which the gas properties vary over the field, was 
similarly first considered by W. Nusselt (ref. 2). The following ex- 
pression is obtained for the heat transfer under the restricting con- 
dition that the temperature dependence on v and a may be repre- 
sented by exponential laws in the absolute temperature, the exponents 
of which for gases of the same substance must be equal^; 

Nu = F(Gr, Pr, Te) (3) 

The form here chosen for the temperature characteristic^ 


^In the notation of the characteristic numbers, the proposals of 
the Heat Conference at Koln in 1931 are followed. See Forschg. Ing.-Wes. 
Bd. 2, p. 380, 1931. 

^ee eq. (5). 

^This is approximately the case, for example, for air, oxygen, and 
hydrogen; see reference 1, p. 433. 

^The abbreviation Te is chosen to conform to the other abbrevia- 
tions for the characteristic numbers and to remind the reader that it 
refers to a temperature relation. 
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Te = 0/T. (4) 

(T^, T«,^(°K) = absolute temperatures of the body and the surroundings, 

respectively) differs inessentlally from the nondlmensional ratio 
T^/Too chosen by Nusselt, namely by the constant -1. Equality of Gr, 

Pr, and Te in two different cases means complete similarity, that is, 
similarity of the velocity field, the temperature field, and the field 
of all the coefficients characteristic of the substance (ref. 1, p. 434) 
and therefore also equality of the Nu nxunber. The possibility still 
remains, however, that similarity may exist at the same time that the 
previously mentioned condition of the power law is not satisfied. 

In the present report the characteristics represented in equa- 
tion (3) will be computed for the heat transfer from heated wires and 
pipes in air, hydrogen, and oxygen; and from these values the most 
probable form of the heat-transfer law, obtainable at present, will be 
determined for the case of the horizontal cylinder in diatomic gases 
(Pr = 0.74).^^ The dependence of Nu on the temperature characteris- 
tic Te, which follows from the similarity theory and is illustrated by 
an example (ref. 1, fig. l) as a third independent variable, will be ver 
if led for several tests in the range of Gr from 10“^ to 10. For the 
smallest values of Gr, the effect of Te considerably outweighs that 
of Gr. In the region of large Gr from 10^ to 10^, on the contrary, 

Te is practically without effect. The large scatter of the test points 
in the nondimens ional representation of the test results of Nusselt 
(ref. 2) and Davis (ref. 3) (±32 percent in the case of Nusselt) is due 
to the fact that the parameter Te is not taken into account in accord- 
ance with equation (l). The dependence of Nu on Gr, Pr, and Te 
will be theoretically clarified heieinafter. The effect of Te for 
small Gr may also be correctly estimated quantitatively. 

From the tests conducted by various investigators for different 
finite Te values, the limiting law of small temperature difference is 
determined in each case by extrapolation to Te = 0. This is of 
primary theoretical and practical significance since it is identical 
with the theoretically simpler case (l) and is thus free from the pre- 
viously mentioned restrictive assiimption with regard to the temperature 
dependence of the constants defining the properties of the substance. 


^®-For the values of X and v required for computing the char- 
acteristics, values are assumed for the temperature t^ (in agreement 

with E. Schmidt). The similarity consideration leaves the choice of 
the reference temperature free. 
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2 . Radiation 

From the measured over-all heat losses directly detenained in the 
tests, the amount due to radiation must he subtracted. For this pur- 
pose, a well known equation (ref. 4, eq. 38, p. 232) was employed; the 
best mean values at present known were substituted as radiation con- 
stants for the metal cylinder surfaces, as follows: for platinum, nickel, 

silver, and tantalum, the values given by H. Schmidt (ref. 5); and for 
copper and iron, the values given by Grdber (ref. 6, p. 196). Small 
deviations as compared with the values used by the investigators them- 
selves are without significance in the range of small values of Gr 
because the radiation component in the case of thin wires amounts to 
only a few percent of the heat of convection. For the tests of Wamsler 
and of Koch on thick pipes, however, for which the radiation losses are 
up to 60 percent, an accurate knowledge of the radiation coefficients 
is required. These coefficients were obtained by Wamsler (ref. 13) 
through his own radiation tests. After correction® of computations of 
reference 13, the values obtained herein for the radiation heat are 
0.2 to 2 percent larger. 


3. Property-Determining Constants of Gases 


In order to compute Nu, Gr, and Pr, values of the following 
properties are required: density p, dynamic viscosity p, heat con- 
ductivity X, and specific heat c for air, hydrogen, and oxygen. 


(a) Constants for air . - The density was computed throughout the 
temperature and pressure range according to the ideal gas law^ . 

The true specific heat for constant pressure was obtained from an 
equation given by Holborn and JakobG for the mean specific heat. 


The temperature dependence of the dynamic viscosity was determined 
by graphical adjustment of the values given by Erk^. These values 


®Wamsler (ref. 13) computes incorrectly with the previously men- 
tioned radiation equation without the factor of the area ratio in the 
denominator. 


^Landolt -Bornst e in : 

5th 

ed. y 

vol. 

I, p. 

43. 

^Landolt -Bornst e in : 

5th 

ed. , 

vol. 

II, p. 

. 1274. 

^Landolt -Bornst ein : 

5th 

ed. y 

vol. 

I, pp. 

. 177-181, Eg. I 


144, Eg. Ila, pp. 138-141. 
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differ only slightly (at most, 4 percent) from the graphical mean 
values given hy Erk from 0° to 700° 

The heat conductivity at high temperatures was obtained "by a gas 
kinetic extrapolation, since the test values collected by Jakob^ lie 
only between 0° and 212° C. The Nusselt formula (ref. 1, p. 491; 
ref. 6, p. 192), which deviates from the available test values by at 
most 2.8 percent, is employed herein. 

1 2 

From the gas kinetic relation-^ 

Pr = p c/X = (n + 2)/(n + 4.5) (5) 

where n is the number of degrees of freedom of the molecular motion 
according to which the Pr number of a gas is independent of the tem- 
perature (and moreover is the same for all gases with the same number 
of atoms); there is obtained for air a small decrease of Pr = 0.739 
for 0° C to 0.721 at 200° C and 0.679 at 1000° C (value obtained from 
formula for n=5 is Pr=0.737), which, however, is of no signifi- 
cance because it lies within the uncertainty limits of the p and X 
values. The same holds for the differences in the Pr numbers among 
air, hydrogen, and oxygen. 

(b) Constants for hydrogen and oxygen . - For the kinematic vis- 
cosity and heat conductivity the values given by Davis (ref. 3) are 
employed, which also for air are in good agreement (deviation, at most 
2.7 percent) with the chosen values. A knowledge of the specific heat 
is unnecessary if the Pr value for 0° C is satisfactory; the values 
are 0.717 for hydrogen and 0.731 for oxygen^^. 

(c) Gas constants for high pressures . - From the tests of Petavel, 
which extend up to 160 atmospheres, only those up to 40 atmospheres 
were computed in order that reliable gas constants would be available. 
In this pressure range and up to 1000° C, the computation was conducted 


^^ien -Harms: Handb. d. Exp. Phys., IV, 4, p. 531. 
^^Landolt-Bbrnstein: 5th ed. , vol. II, p. 1304. 

^^A. Busemann in Wien-Harms: Handb. d. Exp. Phys., TV, 1, p. 359. 
This relation may be obtained from an equation given by Eucken (Physik 
Z., vol. 14 (1913), p. 324. 

^^A. Busemann in Wien-Harms: Handb. d. Exp. Phys., IV, 1, p. 362. 
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with the normal values since the dependence on the pressure of X and 
p and the deviations of p from the ideal gas law are either not known 
at all or are only partly known, and the known deviations amount to only 
a few percent^'^. 


4. Individual Invest igations^^ 

(a) General. - The nondimensional representation of the results is 
effected according to equation (3) hy logarithmic graphs with Nu as 

a function of Gr and with Te as the only parameter, since Pr 
(according to sec. 3a) may be considered as constant. Even if the 
slight variation of Pr there indicated is considered as real, there 
would, according to a test result of Davis (ref. 3) on the effect of 
the Pr number, be obtained only displacements of the curves of the 
order of magnitude of l/lOO of those which can actually be observed. 

In table I are summarized the test objects and test conditions of 
the individual investigators and the range of Te used. In order to 
judge to what extent the condition of infinite extension of the fluid 
was satisfied, the ratio of the height H of the surrounding space to 
the cylinder diameter d is shown. In figures 1 to 7 , approximately 
equal Te values are Indicated by the same symbols. The decrease in 
Nu with increasing Te is particularly evident in the individual 
point groups (because of the same wire), whereas the values for different 
wires naturally show larger scattering as compared with one another. 

The curves Te = constant were obtained by a graphical adjustment pro- 
cess and, similarly, the extrapolation points and curves for Te = 0 
and the dependence of Nu on Te (relative decrease in Nu for 
ATe = l). The indicated minimum and maximum values hold for the indi- 
vidual wires (point groups); the indicated mean values, for the individ- 
ual investigators. 

(b) Results of Ayrton and Kilgour (ref, l), figure 1 . - The authors 
of reference 1 give in the form of diagrams for each wire the total heat- 
transfer coefficient as a function of the wire temperature; two wires 
(0.0206 and 0.0282 cm diam. ) because of the very large scattering and a 
third wire (0.0102 cm diam.) because of its strong deviation in position 
and inclination were excluded from the evaluation of the tests. The 


^‘^Landolt-Bornstein: 5th ed. , Eg. I, p. 64 ff. Eg. Ila, p. 144. 

^^Acknowledgement is made to Dr. K. Winkler, Leipzig, for his 
assistance with the computations. 
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points computed from the low temperatures of 40° to 80° C were simi- 
larly not used on account of the large scattering. 

(c) Results of Langmuir (ref. 8)^ figure 2 . - Langmuir gives in a 
table the adjusted values of the total energy losses as a function of 
the absolute wire temperature increasing in stages of 200^ C. Only 
the tests up to 1300° K were computed herein in order to have reliable 
values available for v and X. The convergence of the curves 
Te = 1.33 and Te = 3.33 shows that the dependence of Nu on Te 
decreases with increasing Gr^ the decrease independently determined 
for the individual wires being: 


log Gr 

0.1-4 

0.8-4 

0.6-3 

0.5-2 

0.4-1 

Decrease, percent 

7.2 

7.5 

5.4 

3.2 

2.2 


This experimental confirmation of the results following from theore- 
tical considerations (sec. 6b) acquires significance in that the tests 
of Langmuir are characterized by a large change of the Te number and 
smallest scattering. The extrapolation from Te = 1 to Te = 0^ on 
account of the absence of test points for this value^ was effected with 
the mean value of 13.2 percent taken from reference 7. 

Langmuir (ref. 9) also discusses his own tests at various air 
pressures (lO to 760 mm Hg) and variable room temperature (t* = -190° C 
to 600° C) which on account of the strong change in the expansion coef- 
ficient (3 = 1/T«,) would be of great interest and the only tests of 
their kind. Unfortunately^ these tests have not as yet been published, 
and the test data are not available in a form in which they may be 
evaluatedlS. 

(d) Results of Bijlevelt (ref. 10), figure 3 . - Since the tests of 
Bijlevelt were primarily for the purpose of investigating forced con- 
vection, the experimental setup was designed for the sensitivity re- 
quirements of this flow. This explains the greater scatter of the test 
points of reference 10 as compared with those of other investigators. 

For each individual wire (with the exception of the Ta and Ni wire), 
a grouping is nevertheless observed with respect to Te in the usual 
sense. 

(e) Results of Kennelly, Wright, and Bijlevelt (ref. 11 ), fig- 
ure 4 . - The investigation of reference 11 is noteworthy in that through 
the pressure (p, v) changes, even for constant wire diameter and constant 
surrounding temperature, a change is effected in the Gr number so that 


^^For these data the author is indebted to a personal communication 
from Dr. I. Langmuir. 
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curves of Te = constant are here obtained from tests on a single 
wire. The nondimens ional representation of these tests (fig. 4) has 
been given previously (ref. l) . The authors give for each wire a dia- 
gram with the heating current as a function of the pressure for three 
or four different temperature differences, as parameter of which the 
values for five different pressures were computed in terms of the 
characteristics and were connected by the curves Te = constant. The 
extrapolation to Te = 0 shows extremely good agreement for both wires. 
The values of the wire of medi\mi diameter (d = 0.02616 cm) were not 
plotted because of the irregularities obtained. 

(f) Results of Petavel (ref. 12), figure 5 . - From the over-all 
heat-transfer coefficients given for constant pressure as a function 
of the temperature difference, values were computed for 200°, 400°, 

600°, 800°, and 1000° C temperature differences up to pressures of 

40 atmospheres. For the sake of clarity, however, only the curves for 
the extreme Te values are plotted in figure 5. Since the heat con- 
ductivity and density of hydrogen differ from the values of air by the 
factors 6 and 14, respectively, these tests would have been particu- 
larly suitable for checking the similarity law. Unfortunately, however, 
these tests are evidently unreliable because of the far too small 
Jacket pipe (H/d = 18.6 as the smallest value of all investigators). 

The similarity law for the three gases is confirmed only for Te = 0.69 
and Gr > lO^. In the case of air and oxygen, there is a splitting 
with respect to the dependence on Te, but in the reverse sense from 
that observed for all wires of the other investigators. In the case 
of hydrogen, a dependence on Te in the correct sense is observable 
only at the smallest values of Gr. Finally, the Nu values also for 
small Gr lie too high as compared with the values of all other in- 
vestigators and the rise in the Nu with increasing Gr has too great 
a lag. For these reasons, the results of Petavel for determining the 
heat-transfer law were not used. 

(g) Results of Wamsler (ref. 15), figure 6 . - Wamsler's curve, 
which was obtained on the basis of correctly computed radiation loss 
(sec. 2), shows a dependence on Te such that the mean values for 

Te = 0.60 lie about 5 percent higher than those for Te = 0.20. This 
small splitting of the effect of Te which occurs in the opposite 
sense to that in the case of small Gr, however, in all probability 
does not contradict the facts but may be explained by the consideration 
that the radiation coefficients in the entire temperature range from 
50° to 270° C were assumed constant, whereas in general for metals there 
has been established an increase of the radiation coefficient with the 
temperature. If in the previously mentioned temperature range a rise 
in the radiation coefficients of about 5 percent is assumed, as follows 
for cast iron from Wamsler's determination of the radiation coefficient 
and for wrought iron from the determination of Nusselt (ref. 14), the 
observed splitting of the effect of Te is neutralized and the curve 
shown, valid for all Te values, is obtained. The test values of the 


CH-2 


MCA TM 1366 


9 


copper pipe and of a wrought iron pipe (of 5.9-cm diam. ) were not em- 
ployed in figure 6 because of the obvious irregularities indicated by 
the scattering, 

(h) Results of Koch (ref. 15), figure 7 . - The nondimensional 
characteristics were redetermined herein with the values X and v 
for wall temperature, since the characteristics given by Koch were com- 
puted for a mean temperature of the substances. The test results show 
an unusually small scattering which is to be ascribed a careful account- 
ing for various factors (e.g., determination of the temperature of the 
air in the room and at the walls, and also the end effect of the pipes). 
For the two intermediate pipes, after graphical adjustment, a decrease 
of Nu for Te = 0.07 to 0.52 by about 2.5 percent is observed, where- 
as in the case of the smallest pipe no ordering is observed and in the 
case of the largest pipe there is a partial decrease and a partial 
scattering. In agreement with the results of Warns ler, the dependence 
on Te is thus practically zero also for large Gr (lO^ to 10^ ) . 


5. Svumnary of Test Results 

The results of the different investigators as regards their agree- 
ment may be compared in two groups. In the region of small Gr 
(10“^ to 10), the tests of Ayrton and Kilgour, Langmuir, Bijlevelt, and 
Kennelly and coworkers are in agreement. From these investigations are 
determined, on the one hand, the curves Te = 0 for the limiting law 
of small temperature difference (fig. 8) and, on the other hand, the 
curves Te = 0.65 corresponding to an intennediate test range (fig. 9). 
For Te = 0 the maximum scatter of the results of the four- investigators 
is 10.5 percent, the minimum 5 percent, the corresponding values for 
Te = 0.65 being 22 and 8.5 percent, respectively. The middle curves 
of figures 8 and 9 obtained for Te = 0 and Te = 0.65 give for this 
change in Te an average change in Nu of 14.7 percent (fig. lO). 

A relative decrease in the Nu number by 22 percent with increase of 
Te from 0 to 1 is thus obtained as the experimental mean value of the 
effect of Te for small Gr. 

In the region of large Gr (lO^ to 10^), the test values of 
Wamsler and Koch from energy measurements and the values obtained from 
the temperature field measurements of Jodlbauer (ref. 16) by integrat- 
ing over the cylinder perimeter are shown in figure 10. As follows 
from the discussion of the results of Wamsler and Koch, the values of 
Koch must be considered as more reliable. There is nevertheless an 
uncertainty in the values due to the uncertainty in the radiation com- 
ponents which lie between 40 and 60 percent. On the other hand, the 
temperature field measurements of Jodlbauer are free from this source 
of error. Since, moreover, the results of reference 16 are in good 
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agreement (mean deviation^ 4 percent) with the theoretical solution 
given in Part II they were assigned a greater weight in the evaluation. 
Jodibauer himself, however, considers his energy measurements, which 
are about 10 percent higher and lie near the values of Koch, to be 
correct and his field measurements in error because of the introduction 
of the thermocouple. The test points of Wamsler are, on the average, 

18 percent higher than the theoretical values, and those of Koch, 

11 percent higher. The interpolation curve is not everywhere satis- 
factory as regards its slope and curvature and could be determined only 
with a partial deviation from the test values. 

The differences still existing, after extensive adjustment, among 
the individual investigators and the difficulty of interpolation show 
that the heat-transfer curve is still not conclusively determined with 
the desirable accuracy of several percent, but that, on the contrary, 
further tests are required. The given interpolation curve nevertheless 
represents the most probable curve of the heat -transfer law, obtainable 
at the present time, for diatomic gases with account taken of the de- 
pendence on Te in the range of Gr from 10“^ to 10^. The numerical 
values of this curve are given in table II. 


6. Qualitative Theoretical Interpretation of Heat-Transfer Law 

(a) Streamlines. - It is necessary to explain theoretically the 
dependence of the Nu number on Gr and Te, determined experimentally 
in sections 4 and 5 (fig. lO), and also the dependence on Pr as ob- 
tained by A. H. Davis (ref. 3) for liquids. A solution of the differ- 
ential equations under the assumption of a thin (as compared with the 
cylinder diam. ) heated laminar layer is given in Part II for diatomic 
gases (Pr = 0.74) and small temperature differences (i.e., approximately 
constant properties, Te » O) and is valid for mean values of Gr of 
about 10^ to 3><10®. On the basis of these theoretical results and with 
the aid of suitable assumptions, the dependence of Nu on small and 
large Gr and also on Pr and Te will be discussed theoretically 
without any necessary implications of finality in the conclusions. 

The streamline plot (fig. 11 ) obtained from the solution and the 
isotherms of figure 12 show that cold air streams from the bottom and 
sides of the cylinder with increasing velocity in laminar flow upwards 
along its surface, is deflected away from the surface in the region of 
the upper stagnation point, and forms a rising current of warm air. 

With increasing Gr of the cylinder, the thickness of the streaming 
layer relative to the cylinder diameter decreases. Since the velocities 
increase at a higher rate, however, the Reynolds mmber Re of the 
boundary layer increases so that the boundary layer finally becomes 
turbulent. According to the schlieren photographs of figure 23, this 
occurs in the neighborhood of the upper stagnation point at Gr = 3.5x10 . 
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With further Increasing Gr, the point of turbulence transition travels 
upstream (fig. 26) and reaches the equatorial region at Gr = 3^10^. 

With decreasing Gr of the cylinder, the Re of the boundary layer 
likewise decreases, and at the same time the thickness of the layer in 
relation to the cylinder diameter increases so that the given solution 
of the differential equations finally (below Gr = 10^) becomes in- 
valid. At very small Gr (very small Re), the isothermal picture of 
the free flow in the neighborhood of the body becomes increasingly 
similar to the concentric isothermal picture of the pure heat con- 
ductionl^ so that almost Sitatic heat- conduct ion relations may here be 
supposed. An estimate of the Te effect made under this assumption 
at small Gr likewise gi^es good agreement with the measurements (see 
following section b). Th5 following fundamental difference is never- 
theless to be observed. |'or initially given cylinder temperature and 
heat transfer, the soluti^ of the heat-conduction equation has as the 
potential function at infinity a singular point (negatively infinite 
temperature). That is, even at a very large distance from the cylinder 
the room temperature continues to decrease, whereas actually in the 
case of free flow the room temperature at large distance from the body 
very soon becomes constant. 

(b) Heat transfer. - Dependence on Gr. For the dependence of Nu 
on Gr, the solution of the differential equations under the familiar 
assumptions of the Prandtl boundary- layer theory gives 

Nu » Gr^/^ (6) 

as is also approximately shown by the tests of Warns ler and Koch (fig. 10). 
At very large Gr with at least partial turbulent flow, the heat trans- 
fer must rise more strongly on account of the increased mixing. With 
increasing turbulent mixing the dependence of the heat-transfer coef- 
ficient on the position must also decrease. The assumption that a is 
independent of the position (or on the cylinder diam. ) , as is obtained 
from the tests of Griffiths and Davis (ref. 17) for a vertical plate, 
gives on the basis of dimensional considerations 

Nu » GrV^ (7 ) 

Above Gr = 3.5X10®, the proportionality of the Nu number to the 
fourth root of Gr corresponding to the upstream travel of the point 
of transition to turbulence must gradually become a proportionality to 
the cube root. Experiments of W. King (ref. 18) on vertical cylinders. 


^^See photographs by R. B. Kennard; Bur. Stand. J. Res., vol. 8 
(1932), p. 787. Report on this by M. Jakob: Forschg. Ing.-Wes., vol. 4 
(1933), p. 45. 
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plates, and blocks at Gr up to 10^^ confinn the cube -root law as 
already pointed out by M. Jakob and W. Linke (ref. 19). For very small 
Gr, for which the amount of heat transferred by convection is always 
below that transferred by conduction, the Nu number in accordance 
with its definition must gradually become constant, as may be seen from 
the discussion of data below Gr = 10“^ (fig. lO) . 

Dependence on Pr. The dependence of Nu on Pr for 
Gr = constant and Te = constant may be understood by assuming a 
change in X with all the remaining magnitudes kept constant. The 
increase in Pr, through a decrease in X, gives in the region of 
predominant heat conduction (very small Gr) a proportional decrease 
in the quantity of heat transferred, that is, Nu = constant. In the 
region of large Gr, a decrease in X (on account of the smaller con- 
duction component as compared with the convection component) gives a 
relatively smaller decrease in a, that is, an increase in Nu. In 
agreement, investigations by Davis (refs. 3 and 20 ) on liquids of very 
different Pr (Pr =0.74 to 7940) for small Gr give the limiting 
value of Nu which is independent of Gr and also of Pr. For medium 
and large values of Gr, an increase of Nu with Pr was obtained. 

Dependence on Te. The dependence of Nu on Te for 
Gr = constant and Pr = constant for diatomic gases is determined by 
the temperature dependence on X, v, a, and also by the temperature 
chosen in computing the characteristics Gr and Nu (Pr is independ- 
ent of the temperature) for the fluid properties X and v (P = 1/Too 
is constant in the entire field, ref. l); the wall temperature t^ or 

T,, was chosen herein, 
w 

An increase of Te = (T^/T^)-1 is considered to be due to a de- 
crease in T with T,, constant so that the values of X and v 

GO W 

employed for the computation of Nu and Gr remain unchanged. The 
associated Increase of 3 and 0 Is compensated hy a decrease In the 
gravity field g so that g30, that Is, the lift acceleration and Gr 
remain constant. Hence colder outer layers with smaller X and v 
hut larger p now take part In the heat transfer. In the region of 
very small Gr with predominating heat conduction X, the magnitude 
a, and hence Nu, must therefore decrease with Increasing Te. This 
Is also shown hy the data discussed previously between Gr = 10“^ and 
10 (see fig. 10 ). For the region of very large Gr with smaller vis- 
cosity and heat -conduction effect as compared with the turbulent momen- 
tum and heat mixing. It also follows^® from the momenta equation. 


^®Gnder the reasonable assumption that with the change of Te the 
profiles of the velocity and temperature to a first approximation under- 
go affine variations. The 'standard velocity W may thus be taken for 
example as the maximum value of the tangential velocity at the equa- 
torial region of the cylinder. 
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because of the constancy of the lift acceleration, that the standard 
velocity W is constant. The participation of colder layers with in- 
creased value of p in the heat transfer here means an increase in the 
heat -transfer coefficient a, which is proportional to cpW, and there- 
fore in Nu. 


The decrease of Wu with increasing Te for very small Gr 
contrasts, therefore, with an increase for very large Gr. The effect 
of Te for small Gr values must therefore, in agreement with the 
test results of Langmuir, first decrease with increasing Gr values 
and finally, somewhere in the range of medium Gr, must practically 
vanish, as is evidently the case according to the experiments "between 


Gr = 10^ and lo"^ (fig. 10 ). 


For the range of very small Gr, the order of magnitude of the 
dependence of Nu on Te under the approximating assiamption of purely 
static heat conduction may he estimated when a mean conductivity is 
used. 


X 

m 




X(T)dT 


( 8 ) 


the following expression is obtained for the Nu niambers corresponding 
to two different Te values Te^^ and Te^ for equal wall temperature 



Nu(Tei) Xj^(Tei) 
NuCTeg) ■ ^mCTeg) 


The value Te^^ =0 is chosen, that is, vanishing temperature difference 
or T^^ = 1^; 8-nd Te2 = 1, that is, T^^/T^^ = 2. For the temperature 

dependence of X, the previously (ref. 1) given exponential form is 
usedl9 


with n = 0.738^° 


X/Xq = (t/Tq)^ 


( 10 ) 


To- 



is the value of 


X 


for the arbitrarily chosen temperature 


Of) 

^Determined from the previously given exponent for a of 1.738 
neglecting the temperature dependence on c. 
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From 


^m(Te2) 



X dT 


with Tq = 


in equation (lO) 


and hy integration: 


\n(Te2) 


>^(Twi) 
n + 1 




1 


Since Xjj^(Te3^) = X(T^^), the following expression results after sub- 
stitution of the numerical values: 


2 f, /ly-'^SS 
X^(Te^) - 1.738 | ~ \2J 


0.81 


( 11 ) 


The decrease of Nu for Gr = constant with increase from Te = 0 
to Te = 1 is thus theoretically obtained as 19 percent. The pre- 
viously discussed tests of the four investigators between Gr = 10-^ 
and 10 gave values between 13 percent and 26 percent with a mean 
value of 22 percent (see fig. 5). This good agreement justifies the 
assumption of approximately static conducting conditions at small Gr. 


II. THEORETICAL SOLUTION OF THE BOUNDARY EQUATIONS FOR THE 
HORIZONTAL CYLINDER (CASE OF STEADY MOTION) 

1. Abstract 

TLe large number of experimental investigations of heat transfer 
in free convection correspond to only a single physically satisfactory 
theoretical solution of the differential equations, namely the solution 
for the vertical plate (two-dimensional steady-flow case) for moderate 
temperature differences, which was given by Schmidt and Beckmann 
(ref. 21) with the help of Pohlhausen. The solution is based on the 
approximation of the Prandtl boundary -layer theory (laminar flow in a 
layer which is thin as compared with the distance from the lower edge; 
velocity normal to the wall small as compared with that in the direction 
of the principal flow along the wall) . The results of the theory as 
regards the velocity and temperature fields and therefore the heat 
transfer agree very well with the corresponding measurements. 
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For the case of the horizontal cylinder, schlieren photographs in 
air (see ref. 22 and figs. 23 and 25) show that, for Grashof character- 
istics (Gr = d^gP0/v2) of about 10^ and above, the heating of the air 
extends to only a relatively thin layer around the cylinder, a fact 
which permits the mathematical simplification of the problem. On the 
assumption that the heat-transfer and the flow processes are restricted 
to a thin film, as compared with the cylinder diameter, with laminar 
flow (boundary-layer ass\xmption) the differential equations (ref. 23) 
for moderate temperature differences can be solved approximately for 
the velocity and temperatiire fields; and therefore the heat transfer 
can be computed which, in the range of Gr 10^ to 3KL0°, within which 
the Initial assumptions are satisfied, is in good agreement with 
experiment . 


2. Setting up of Differential Equations 


To start, the hydrodynamic differential equations are used in a 
form previously given (ref. l)^^ which shows a "lift term", charac- 
terizing the free convection, that arises from the combining of the 
gravity term and the hydrostatic component of the pressure drop, so 
that only the gradient of the dynamic pressure p* remains as the 
pressure force. This dynamic pressure is equal to zero where there 
is no motion and no temperature difference with respect to the medium 
at a large distance. As was indicated in detail in reference 1 (p. 
429), the density ant^ other characteristics of the substance in the 
entire field may be considered as constant (case of small temperature 
difference) snH the equations will continue to describe the free- 
convection problem* The equations in vector form are then^ 


wO grad w = - l/ p grad p* - v rot rot w 

div w = 0 
grad 0 = a A 0 


- gp0 




) ( 12 ) 


y 


with the boundary conditions: w = 0, 0 = 0 on the cylinder surface; and 

w = 0, 0 = 0, p* = 0 at infinity. (w, velocity; 0, temperature differ- 
ence with respect to temperature at infinity; p*’ , dynamic pressure; 
p, density; v, kinematic viscosity; 6, expansion coefficient; a, temper- 
ature conductivity; g, gravity acceleration vector; g, its absolute 
value; r, cylinder radius.) 


^^In equation (3c) in reference 1 a minus sign mistakenly appears 
on the right side. 

^^In place of Aw the invariant form is written grad div w - 
rot rot w. 
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For treating the two-dimensional flow about the cylinder, these 
equations are conveniently written in arc length coordinates with the 
arc length s on the cylinder perimeter (taken in clockwise direction 
with the lower stagnation point as origin) and the normal distance n 
from the cylinder surface, and the velocity components u and v 
par-allel and normal to the wall, respectively. There is then obtained, 
without neglecting any terms^^: 



du 


5u 

— u 

- n 

SF 

+ V 


1 

• 1 

e 

•V' 

d^u 


1 

5u 

dn2 

4- 

r 

+ n 

5n 

l/e2 


1/e 


r 


Sv 

1 \r 

r + 

u 

n 

5s 

~r V 


r + n 


1 Sp* r • s 

— + gB0 sin — + 

p ds r + n r 


d"=u 


u 


2r dv 


u‘- 


§n r + n 
1 • e e -1 1 


1 dp* 
p "ST 


- gpe cos - + 


jd^v 1 dv 

V < — ^ + + 


d2\ 


2r du 


dn^ r + n dn (r + n)^ ^^2 (j. + n)‘ 


6^ i/e 


(r + n)^ ds 
2 • 1 


r du dv V _ 

+ n d¥ Sn r + n 

lie 


(13) 


u 


de , 

de J 

d^0 1 1 d0 1 r2 


Si 

dn 

. 2 r + n (j. + n)2 

ds2 

1 

e-l/e e2 

'1/62 i/e 

1 


The approximations made according to the boundary -layer theory are now 
introduced. The cylinder radius r and hence also s (with the ex- 
ception of the lower stagnation point) and u are assumed to be of the 
order of magnitude 1; n is assumed to be of the order of magnitude of 


^^he hydrodynamic equations agree, except for the dynamic pres- 
sure and lift terms, with the equations (lO), (ll), and (l2) for 
r = constant of W. Tollmien: in Wien-Harms, Handb. d. Exp. Phys., IV, 
1, p. 248. 
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the "boundary- layer thickness £ so that n + r may be re- 

placed by r. The estimate of the order of magnitude of the individual 
terms according to the method given by Prandtl and his coworkers then 
gives the results written underneath the individual terms. In order 
that the friction and inertia forces may be of equal effect, it is 
necessary that v (and correspondingly a) be of the order of magnitude 
of e^. In order that the lift forces for the motion in the direction 
of the main stream be of significance as compared with the inertia and 
friction forces, that is, in order that 'free convection' exist at all, 
it is necessary that gP0 » 1. From the two momentum equations it 
then follows that the pressure terms are at most » 1. Since p"*" is 
zero at the outer edge of the boundary layer, it follows from 


by integration from the outer edge of the boundary layer to the cylinder 
surface over the length e that ^ /p in the boundary layer can at 
most be of the order of magnitude e, that is, so that l/p dp*/Ss = e. 

Since in what follows only the terms of the order 1 are retained, 
the following result, important for simplifying the mathematical treat- 
ment, is obtained: The tangential drop of the dynamic pressure is to 

be neglected as compared with the tangential lift, inertia, and friction 
forces. The term p* alone is contained in the equation of the normal 
momentum. Hence, u, v, and 6 are to be computed from the equation of 
the tangential momentum, the continuity equation, and the heat-transfer 
equation (eq. (l4)). Thereafter p* may, if desired, be computed from 
equation (15) of the normal momentum. Hence, the boundary -layer equa- 
tions for the horizontal cylinder for free convection are 


1 hv* 

p dn 



\ 



( 1 ^) 


1 Sp* nn 

± gP0 cos 

r p on 


r 


s 


( 15 ) 
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A corresponding estimate of equations (13), first rendered non- 
dimensional, gives the two corresponding conditions: 

_ . and 1, i.e. Or - ^ (b) 

For a maximum e, in the sense of the approximations, of 0.1 (approxi- 
mation to about 10 percent), this gives Re = 100 and Gr = 10^ as 
lower limits of the computations. The approximation is therefore better 
the larger Gr, that is, the smaller the boundary-layer thickness e . 

The upper limit is then attained if the boundary layer becomes turbu- 
lent, which (according to part IV) is the case for Gr = 3.5><10® at 
the upper stagnation point, for Gr = SxlO^ at the equator. The range 
of applicability of the laminar -boundary-layer computations therefore 
lies between Gr = 10'^ and Gr = 10^. 


3. Reduction to Ordinary Differential Equations 


The usual introduction of the stream function 


u 






(16) 


eliminates the continuity equation. By introducing nondimens ional 
magnitudes with Gr' = r3gp0/v2 -by means of the transformations^ : 

X = ^ y “ ? Gr'lA Y = ^ Gn'-lA ^ g 


there remain of the five constants (v, gP, a; 0, r) of differential 
equations (14) only two in the nondimensional combination of the Prandtl 
characteristic Pr =v/a. 


ay 

Sy ax ay " 3x ^y2 ^y3 


T sin X 


(18) 

^ _ i_ 

Sy Sx Sx by Pr Sy2 


Equation (l7) is not an arbitrary transformation which carries 
equation (14) over into equation (l8)^ but a definitely determined 
transformation with the aid of which it is possible to eliminate Gr 


^"^Gr* referred to r is more convenient for the computation than 
Gr referred to d. In the final result^ the more usual Gr is again 
introduced. 
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so that it no longer occurs in equation (18). The solution of equa- 
tion (18) therefore likewise no longer contains Gr. It is therefore 
unnecessary to know the solution of equation (18) if it is of interest 
to know only how the solution (e.g., velocity, temperature, houndary- 
layer thickness, heat-transfer coefficient) of the initial equation (14) 
depends on the Gr niomber. This dependence of Gr is already com- 
pletely represented by the transformation (17), whereas equation (18) 
contains the further dependence on Pr and on the space coordinates. 

In other words, the solution of equation (18), obtained for a definite 
Pr value24a gives, by means of equation (17), at the same time the 
velocity and temperature fields for all Gr values. On the other hand, 
the transformation (17) without the solution of equation (18) already 
gives important general information on the flow condition and the heat 
transfer for free convection with respect to the dependence on Gr. 

From equation (l7) there follows, if the following expression is set up. 


U 


St 



(19) 


that 


r. t 

- = y Gr' 


-1/4 HI = 

V 


U Gr 


'1/2 H^vGr'1/2 l|®=|lGr'l/^ 

V 0 dn dy 

( 20 ) 

That is, the relative boundary -layer thicknesses .decrease with Gr' -1/4. 
the nondimens ional velocities increase with GrV^; and the non- 
dimensional temperature drop increases with Grl/^, from which there 
follows directly the well-known l/4-power law of heat transfer in free 
convection: 


Nu » GrV^ (21) 

These theoretical results agree with those for the rectangular 
plate and are not restricted to the latter and the horizontal cylinder 
alone, but depend on more general assumptions (ref. 19). They are 
valid, in general, for all two-dimensional cases of free flow with a 
boundary-layer character (ref. 23), that is, where the thickness of the 
heated laminar layer is small compared with the distance from the inci- 
dence edge and the radius of curvature of the wall, and the latter does 
not change dis continuously in the flow direction. The existence of 
sharp edges about which the flow occurs is therefore excluded. In this 
more general case, there are on the right side of equation (13) addi- 
tional terms which are determined by the radius of curvature and which 


2^^It is here assumed as self evident from physical considerations 
that one and one only such solution exists. On the mathematical side 
of these existence and uniqueness proofs, see the discussion in refer- 
ence 1 (p. 428). 
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vary with the arc length (dr/ds) hut which drop out in the estimate 
made under the initial assumptions. Furthermore, in place of the sine 
and cosine of the lift term there appear functions which determine the 
direction of the surface element with respect to the horizontal and 
vertical hut which similarly do not change anything in the transfor- 
mations (17) and the conclusions from equations (20) and (2l) based on 
these transformations. 

In agreement with these theoretical conclusions equation (21 ) has 
alreadv been confirmed for vertical plates (ref. 21), horizontal cyl- 
inders^S^ and cylindrical layers (ref. 24), and even for the heat 
transfer from the two sides of a square horizontal thick plate with 
edges (ref. 25) and for the heat transfer with evaporation at a verti- 
cal cylinder in water and carbon tetrachloride (ref. 26) for which the 
previously mentioned assumptions are not all satisfied. 

A reduction of the partial differential system (18) to systems of 
ordinary differential equations is attempted by setting 

q = y g(x) (22a) 

Y(x,y) = p(q)f(x) (22b) 

T^(x;,y) = t(q) (22c) 


for example^ 

U = g| = f(x)g(x)p'(q) 


(22d) 


based on the notion that all the profiles of the stream function f(y), 
tangential velocity U(y), and temperature x(y) for the different 
azimuths x of the cylinder are produced by affine distortions by 
means of the azimuth functions f(x) and g(x) from a single "base 
profile" of the stream function p(q) or the velocity p' (q) and one 
base profile of the temperature t(q) (q = normal coordinate of the 
base profile). Substitution of equation (22) in equation (18) gives 
(primes of p and t denote differentiation with respect to q; primes 
of f and g denote differentiation with respect to x) : 


p'2(f2ggt + ff>g2) - pp"ff'g2 = p'"fg3 + t sin x 
l/Pr gt" + f'pt' = 0 


(23) 


25 


'See Part I, section 5, figure 10. 
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A separation into functions of the two independent variables x 
and q, that is, a successful application of the assumed expression 
(eq. (22)) for solving the equations, is possible only if the two func- 
tions f (x) and g(x) can be uniquely determined from the four ordinary 
differential equations 


f'(x) = ag(x) 

(24a) 

f (x)f ' (x)g2(x) = c sin x 

(24b) 

f^(x)g(x)g' (x) = b sin x 

(24c) 

f(x)g3(x) = d sin x 

(24d) 


(where a,b,c,d are initially undetermined constants) in spite of the 
redundancy of the equations^G. The possibility of such unique determi- 
nation, which is the necessary assumption for the applicability of 
equation (22), will be discussed in more detail after equation (29) is 
considered. If such determination is possible, there remain for p(q) 
and t(q) two ordinary differential equations: 

(b + c)p'2(q) - cp(q)p"(q) = dp'" (q) + t(q) 
t"(q) + a Pr p(q)t' (q) = 0 

In section 4, equation (24) will be solved and in section 
will *be discussed. 


(25) 

5 equation (25) 


4. Determination of the Azimuth Functions F(x) and G(x) 


(a) Setting up of the two differential equations for F(x) . - The 
boundary values of f(x) and g(x) which are required for solving equa- 
tion (24) are determined from physical considerations. From 


St 

= 


dt(q) 

dq 


g(x) 


it follows that 


g(x =0) = go t 0 


(26a) 


^^The author is indebted to Dr. A. Naumann, Leipzig, for several 
mathematical suggestions. 


22 


NACA TM 1366 


Otherwise, the temperature of the cylinder surface would also he the 
temperature in the field along the entire normal at the lower stagna- 
tion point, which cannot he the case. From equation (24a) there is 
then obtained 


f»(x = 0) = agQ (26h) 

Since the tangential velocity at the lower stagnation point must, for 
reasons of symmetry, he different, there follows from equation (22d): 

f(x = 0) = 0 (26c) 

The system of equation (24) with the boundary conditions (26a), (26h), 
and (26c) thus contains five free available constants a,h,c,d, and gQ, 

which are determined in subsequent computation. 

By introducing the normed functions^^ : 

F(x) - G(x) = iM (27) 

ago Sq 

equation (24) may be transformed to: 

F’ = G (28a) 

G G' = ^ sin x (28b) 

F F» g 2 = ^ sin x (28c) 


F g 3 = — ^ sin X (28d) 

ago 


with the boundary conditions 

F(x = 0) = 0 F’(x = 0) = 1 G(x = 0) = 1 (29) 

A necessary assumption for the applicability of the expressions of 
equation (22) for solving the equations is that the two functions F 
and G can be uniquely determined from the four equations (28a) to 
(28d). This would be the case if it could be proven that the four 


27 

By normed is meant fully determined also with respect to a 

factor. 
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equations are not independent of one another^ tut that two equations 
follow from the other two. This proof is partly possible only inso- 
far as one equation (for suitable choice of the constants) is a con- 
sequence of the remaining equations. Thus, elimination of G with 
the aid of equation (28a) gives 

pt pTT _ ■ sin X (30a) 

^0 


F F^^ = 


a^g ^ 


sin X 


F F* ^ 



sin X 


(30b) 

(30c) 


Equation (30c) is thus a consequence of equation (30b) in the case that 
the following first condition for the constants is satisfied: 


d = c/a 

Setting for briefness 

(31) 

c/a2g^‘i = a 

(32a) 


(32b) 

there then remain for the determination of F 

the two equations 

F F* 3 = a sin x 

(331) 

F^ F* F'' = P sin x 

(3311) 


with freely available values of a and p. As a necessary assumption 
for the applicability of equation (22), there remains the further con- 
dition that the two equations I and II possess two equal (or at least 
approximately equal) solutions. 

This proof will be obtained by solving equations I and II by a 
power-series development above and below the equator x = j(/2. For 
this purpose, with (^ = x - jt/2, I and II are transformed into 

F(cp)F^^((p) = a cos (p (341) 

f2((p)F’ ((p)F"(<p) = p cos (|> (3411) 

where the primes now denote differentiation with respect to 9 . 
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(b) Solution of differential equation (l) for F . - The value of 
a, by the boundary values of equation (29), is first determined as 

a = 1 (32c) 


as is seen from the first term of the power-series development of 
equation (33l) about the point x = 0. From the expression for solving 
equation (34l); 

F^(<>) = aQ(l + + Cg <^2 + . . .) (35) 


there follows 

pl((^ = 0) = Bq F^' (cp = O) = aQ C]_ (36a) 

and by substituting in equation (34l) 

aQ^ = a = 1 (36b) 


Substitution of expression (35) in differential equation (34l) 
gives the unknown coefficients, which are taken up to the fifth, in- 
clusive, and which can be successively determined in terms of c-]_: 


Co = 


1 

6 


- ^1^ 


C-t: = - 


18 


54 


+ 7=-T 


"4 = 54 


72 1 1 


^5 - “ 360 '^1 ■ 324 *^1^ ■^216 ‘^l” 


13 

216 


y(37) 


The boundary condition F^(q> = -rt/2) = 0 according to equation (29) 
gives for the value of c^_ the solution of an algebraic equation of 

the fifth degree which is solved by trial 


cj_ = 0.581 


(38) 


Of the five possible real roots it is necessary to choose that one for 
which no further zero of F^ is obtained between cp = 0 and <p = -jt/2. 
The coefficients then assume, according to equation (37), the values 

cn = + 0.581 ^2 = - 0.05626 Cj = - 0.01412 "j 

(39) 

C4 = - 0.00165 C5 = - 0.00066 J 

According to equations (36a), and (36b), there is obtained 

fI(cP = 0) = aQ = 1.504 F^'((p= O) = aQ c^ = 0.874 (40) 
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and the numerical solution of F^(<p) is thereby completely given (see 
table III) • 

It may be remarked^^®' that for equation (33l) a solution was also 
obtained by a power-series expansion at the point x = 0 in powers of 
sin (x/2) which, because of the fact that the curvatiare of F at that 
point is zero, converges with particular rapidity: 

F*(x) = 2 sin I (l + ^ sin2 | + ^ sin4 | (41) 

This solution assumes at the equator the values 

F*(x = jt/2) = 1-500 F*'(x = rt/2) = 0.874 (e) 

which are in good agreement with the values of equation (40 ) of the 
series expansion at the equator. In this way, the sufficient conver- 
gence of f 1(<P) with the five computed terms is assured, and a good 
check is obtained against computation errors. The solution F* will 
not be used in what follows; this is because of the singular behavior 
of equation (33Il) at the point x = 0, and the identity proof must 
therefore be carried out by the series expansion at the equator. 

(c) Solution of differential equation II for F . - For equa- 
tion (3411), corresponding to equation (35), the following expression 
is assumed: 

f11((p) = bo(l + di<P + d2<P^ + . . .) (42) 

Substitution in the differential equation gives the following values 
for the constants: bQ and d^_ are arbitrary 


2bo^d^ 

d^ dj = - I (<li^<l2 + <^2^) 

where, on account of the length of the expressions, d^ and d^ are 
not given as functions of d^^ and d£. The function thus depends 



^^^'Acknowledgment is made to Prof. W. Muller, Aachen, for the 
suggestion that the solution of equation (331) may be reduced to two 
quadratures by separation of the variables. One of the quadratures is 
not in closed form, however, but can be carried out only by a series 
expansion so that this method would not be essentially different from 
the one employed. 
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on the as yet freely available values of Tdq, dj, and 6. The solu- 
tions and are equal if the corresponding coefficients of the 

series are equal. Therefore, the as yet undetermined coefficients are 
equated as follows: 


ho = Eq 

The condition dg = Cg gives 

3 


di = Cj_ 




2bo^di 6 

or with equations (36h) and (44) 

P = - l/3 a = - 1/3 

All the coefficients may now be expressed in terms of d^j^: 


<12 - - I ^1^ 


^5 = 


d^ + 72 di - 72 d;L 


(44) 


(f) 


(45) 


^5 ■ 540 '^1^ 216 

Comparison of the coefficients of the two series for and F^^ 

gives 

bo - ao = 0 d]^ - C3_ = 0 dg - Cg = 0 d3 - C3 = ^ c^^ 


(46) 


d/1 - C/ — - 


216 1 


^5 ■ + 360 ‘^l 405 ^1^ 


(47) 


This result states that the solutions of the two different differen- 
tial equations (34l) and (3411) agree at the equator in the value of the 
function, in the tangent, and in the curvature, hut not in the higher 
derivatives in which they differ fundamentally to a small degree. 
and are approximately, hut not rigorously, equal. The system (24) 

is, strictly speaking, overdetermined and has only approximately unique 
solution. This solution of expression (22) is not rigorously, hut 
only approximately, valid. Because of the agreement of F-^ and F^^ 
at the equator up to and including the curvature, however, these differ- 
ences become more appreciable in the neighborhood of the lower and 
upper stagnation points where the initially assumed houndary-layer 
assumptions are not satisfied. The method that has here been employed 
thus gives an approximate solution of the houndary-layer equations with- 
in the accuracy limits determined by the physical approximations of the 
houndary-layer assumptions. 
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The numerical computation of equations (44) and (46) with aQ 
and cp according to equations (40) and (39) gives: 

bQ = 1.504 d^ = 0.581 ^ " 0.05626 <13 = + 0.01816 


Table III shows the values of the functions fI(<>) and FH((p) computed 
according to equations (39), (40), and (48), with sufficient agreement 
between x = 30°, 150° and even at 165° (graphically interpolated for 
F^^). The negative value of F^^ for x = 0 has no significance 
since, on account of the singular behavior of equation (3411) for 
F = 0, the series (42) beyond this singular point no longer represents 
the solution of the differential equation. For the further equations, 
F^ is employed because only this function possesses the required 
property of vanishing at the lower stagnation point. 

(d) Detemiination of G . - The term G is determined according 
to equation (28d) where, according to equation (49), d = agQ'^. 

Table III shows the computed values based on F^. 


The constants b, c, d which occur in equation (25) can now, with 
the aid of equations (3l), (32a), to (32c), and (45) be expressed in 
terms of a and gQi 


= -0.00321 d 5 = + 0.00144 


(48) 


5. Determination of Basic Profiles p(q) and t(q) 


b = - 1/3 a2 £q4 c = a2 gQ4 d = a gQ4 


(49) 


Equation (25) then becomes 



(50) 


t" + a Pr p t' = 0 


with the boundary conditions: 


q = 0 :p =0 p '=0 t = l 



t = 0 


(g) 


Equation (50) for a = 1 and gQ = 1 may now be numerically integrated. 

The following values were chosen for the as yet free remaining constants 
a and gQ 
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a = 3 gQ = 


(51) 


equation (50) becomes 


p"' + 3pp" - 2p'2 + t= 0 
t" + 3 Pr p t' =0 


(52) 


This system with the same boundary conditions and for Pr = 0.733 
(diatomic gases) has already previously been numerically integrated by 
E. Pohlhausen for the vertical plate so that, because of the introduc- 
tion of the normed functions F and G and the consequent availability 
of the constants, a new numerical integration is unnecessary. The 
stream function p(q) and temperature function t(q) correspond to the 
functions ^(C) and d(^) of Pohlhausen. Their values are given in 
reference 21 in tabular form (ref. 4, p. 187). Graphical representa- 
tions of the velocity profile P* (q) and of the temperature profile 
t(q) are given in figures 13 to 16 (continuous curves). 


6. Complete Solution for Velocity and Temperature Fields 

With the values of equations (5l) for a and gQ, the azimuth 
functions f(x) and g(x) according to equation (27) are now uniquely 
determined. For 

a gQ = 3^4 = 2.280 gQ = 3"V^ = 0.760 (53) 

there is obtained 

f(x) = 2.280 F(x) g(x) = 0.760 G(x) (54) 

Table IV shows the values of f(x), g(x), and f(x)-g(x) computed from 
F^ and G of table III. Figure 17 shows these values graphically. 

With the values f(x) and g(x) from table IV, the values p(q) 
and t(q) from the table of Pohlhausen and the transformations (16), 
(17), and (22) there is obtained the complete solution of the boundary- 
layer equations (14) for the velocity (u and v) and temperature 
fields (0) obtained without any empirical value. Collecting results 
yields (see also eq. (l9) and (20) for simultaneous passage from Gr' 
to Gr ( Gr = 8Gr ' ) : 
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8 



(55a- 
) 55d) 


with X = — and q = 
r ^ 



For the detemination of the normal velocity v it is required to 
know f ’ (x) and g’(x). The term f*(x) is^ according to equation (24)^ 
determined hy means of g(x); (x) was determined hy graphical differ- 

entiation of g(x)^^. Here the qualitative description of the nomal 
velocity obtained from the streamline picture is sufficient. It can he 
shown that the boundary conditions of u^ v^ 6 for q = 0 are cor- 
rectly assumed. If the assumption n <kt which was made with respect 
to the solution of equation (14) is applied for very large q^ there 
are obtained for u and 6 the correct value zero; but for v, in 
contradiction to the boundary conditions (l2)^ a finite value is ob- 
tained which is related to the fact that the cylinder curvature was 
neglected. 

The boundary -layer thickness^ tangential velocity^ temperature, 
and heat-transfer coefficients according to equation (55) behave as 
follows. With regard to the dependence of these magnitudes on Gr, 
there hold first of all the considerations of section 3 in connection 
with equations (20) and (2l). The dependence on the cylinder azimuth 
is such that all profiles of the tangential velocity and temperature 
are obtained by affine distortions from the basic profiles p’ (q) and 


The value l/g(x) (fig- 17) represents the ’extension’ of the tem- 
perature and velocity profiles normal to the surface with increasing 
azimuth, that is, the development of any characteristic distance from 
the wall, for example, of the place of maximum velocity (q = 0.95) or 
of the boundary-layer thickness to be defined later (q = 2.18). In 
accordance with this, the boundary -layer thickness at the lower stag- 
nation point possesses a finite approximately constant value over a 


t(q) . 


values are not given here. 
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large azimuth range hut increases at first gradually then more and more 
rapidly and attains at the upper stagnation point a theoretically in- 
finitely large value (upward current of warm air). The function g(x), 
since it is proportional to the temperature drop, likewise represents 
the variation of the local heat-transfer coefficient with the cylinder 
azimuth. At the lower stagnation point, g(x) is constant and with 
constant slope connects with the values of the other side of the cyl- 
inder. At the upper stagnation point (x = it) this is no longer the 
case; g(x) has a sharp peak which is associated with the fact that the 
houndary-layer assumptions at this point strongly deviate from the 
actual conditions (normal velocity large as compared with the tangen- 
tial velocity) . Under actual conditions this peak is balanced out hy 
a very sharp minimum value as is shown hy the value of Nu along the 
cylinder perimeter computed from the schlieren pictures of E. Schmidt 
(ref. 22)^^. No quantitative comparison of this experimental curve 
of g(x) with the theoretical curve given herein can he given because 
the conditions of the two-dimensional problem are so little satisfied 
for the schlieren pictures (pipe length = 2 times cylinder diam. ) that 
the computed heat-transfer coefficient is 40 percent greater than that 
of the tests of Koch (ref. 22). 

The development of the tangential velocity with the azimuth, that 
is, the increase of any characterizing velocity, for example the max- 
imum velocity, is described by f(x)- g(x) (see fig. 17). The velocity 
increases from the value zero at the lower stagnation point (up to 60° 
approximately linearly) to a maximum at about 128° and, on approaching 
the upper stagnation point (convergence of the flow on either side) 
rapidly decreases to zero. This physically correct behavior of the 
solution at the upper stagnation point is not introduced as a boundary 
condition but is obtained as a necessary consequence of the theory. 

An over-all picture of the velocity and temperature fields is 
given by the streamlines and isotherms (figs. 11 and 
computed from the equations 

= f(x)Gr-l/4 p [E Gr'lA g(x)] = 

= t Gr-lA g(x)] = const 
for various values of the constant s^O, The scale of 


12) . They are 


const 


(55e,f) 


the representation 


^^See figure 13 for maximum Gr = 16x10®. 

®®The streamline picture for azimuth intervals of 15° and in addi- 
tion for 5°, 10°, 170°, 175°; the isotherm picture for azimuth intervals 
of 30° and in addition for 165°, 170°, 175°. 
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holds for Gr = 10^; for Gr = 10® the distances from the cylinder 
surface should he reduced hy a factor of 3.16 and for Gr = 10®, by 10. 
At the lower and equatorial neighborhood of the cylinder up to an azi- 
muth of about 105^, the flow near the wall is directed toward the sur- 
face, while above it is directed away from the surface upwards. In 
this respect the theoretical solution for the cylinder differs funda- 
mentally from that for the vertical plate, which gives only a flow 
toward the plate. The isotherms at the lower part form almost con- 
centric circles. At the upper stagnation point (upward warm current) 
the isotherms theoretically go toward infinity, whereas actually (corre- 
sponding to the gradual dissolution and spreading out of the warm up- 
streaming air) they close at a relatively large distance over the 
cylinder. The isotherm 0.53 gives at the same time the place of the 
maximum tangential velocity. 


7. Heat-Transfer Law 

From the expression for the quantity of heat dQ passing from 
an element of area df, 

dQ = X ( - ^ [n = 0]j df (56) 

the value of the temperature gradient at the wall obtained from 
equation (55c) 

II (n = 0) = ® g(x)t'o (57) 

(t'Q = t'(q = O) = -0.508), and the defining equation for locally 
variable heat-transfer coefficient a(x) 

a(x) = i ^ (58) 

there is obtained 

‘ ^ = g(x)(-tQ*) Gr'V^ (59) 

From this there is obtained for the dependence of the nondimens ional 
local heat-transfer coefficient on the cylinder azimuth x and the 
Gr ntimber (Gr^ = 8Gr')^^^: 


^^Since several formulas of this section will be used later in 
Part III, and in order to avoid misunderstandings, indices (e.g., d,h,H) 
for the characteristics Gr and Nu are written which indicate the 
lengths with which these magnitudes are formed. 
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‘ - = NU(3 ^(x) = 0.604 g(x) 


(60) 


see figure 18. 


To compute the mean heat-transfer coefficient averaged over the 
cylinder perimeter^ 




dx 


there is required the mean value 


(61) 



I g(x) dx; 

Jo 


(62) 


which, by planimetering g(x) (fig. 17), is found to be 

g = 0.616 (63) 

It may also be computed by the following relation from equation (24a) 
with a = 3 according to equation (51) 

3 rt g = f(n) (64) 

which gives g = 0.620. However, the value 0.616 is preferable because 
this value takes into account the total curve g(x) and does not depend 
on the end point of the series expansion f(x). From equation (60) the 
following is then obtained for the mean heat-transfer coefficient as a 
function of Gr for free flow at a horizontal cylinder 

= NU(i = 0.372 Gr^V^ (65) 

A. 

A comparison of this theoretical heat-transfer coefficient for diatomic 
gases (Pr = 0.74) with test results is given in Part I, section 5; see 
also figure 10. 


8. Comparison of Theoretical Velocity and Temperature 

Fields with Available Measurements 

The comparison of the computed velocity and temperature fields is 
carried out with the measurements of Jodlbauer (ref. 16) which are the 


CH-5 


NACA TM 1366 


33 


only ones thus far available^^. Of the fields measured for six differ- 
ent Gr values, the following ones were computed; 

2r = 9 cm, t^ = 99.2°, t = 18.1*^, Gr = 3.76x10® (figs. 13, 14) 

2r = 5 cm, t^ = 104.6°, t = 18.1°, Gr = 6.54xio5 (figs. 15, 16) 

The best agreement with theory is to be expected for the maximum Gr. 

On the basis of the theory it would be desirable to have measurements 
for the largest possible laminar Gr, that is (see Part IV) values of 
about 3.5x10®, which could be realized in air at 20° C at a surface 
temperature of 100° C with cylinders of 42 centimeters diameter. 

The comparison with the theory was carried out in such manner that 
all velocity profiles u(s,n) and temperature profiles 6(s,n) corre- 
sponding to a value of Gr were, with the aid of equation (55) , re- 
computed for the theoretical basic profiles p'(q) and t(q) in the 
representation of which they must all coincide. 

In the case of the velocity profiles (figs. 13 and 14), there is 
a regular deviation in that the measured velocities are greater than 
the computed ones, the deviation from the theory in the maximum velocity 
amounting on the average to 22 percent for the smaller and 17 percent 
for the larger value of Gr. A part of the deviations occurring for 
large distances from the wall beyond the maximum value should be 
ascribed to the uncertainties in the difficult measurements of such 
small air velocities, as can be seen from the following discrepancies. 
The measured velocities at the azimuth 60° for Gr = 3.76x10® lie 
considerably above the theoretical curve, whereas for Gr = 6.54x10® 
they lie on the theoretical curve. On the other hand, the velocities 
measured at x = 30° for Gr = 6.54x10® show strong deviations upward, 
whereas for Gr = 3.76x10® they are in agreement with the theory. The 
measured distance from the wall of the maximum value of the velocity is 
everywhere in very good agreement vfith the theory. 

The measured temperature profiles show very good agreement with 
the theory even for the azimuth of 165°. A small regular deviation is 
noted for all azimuths and Gr numbers for medixim distances from the 
wall. The small deviations at small distance from the wall can be 
ascribed partly to the fact that the condition for the isothermal sur- 
face was not quite satisfied in the tests. The measured temperature 


^^Acknowledgment is made to Prof. E. Schmidt and Dr. K. Jodelbauer 
for their aid in providing the data. 
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drops at the wall for Gr = 6 . 54 x 10 ^ are on the average 3 percent 
smaller than the values given by the theory, while for Gr = 3.76 xL0° 
they are 1 percent higher^^. 

It is likewise of interest to compare the solution based on the 
boundary -layer theory for the vertical plate with the corresponding 
test res\ilts. Reference 21 (figs. 20 to 23) also shows the velocity 
fields to have stronger deviations than the temperature fields. For 
the small plate with smaller deviations, the measured velocities also 
lie somewhat higher than the theoretical values; but for the larger 
plate with larger deviations the contrary is true. For the temperature 
fields, particularly for the large plate, no regular deviation can be 
established between theory and experiment. The deviations occurring 
for the plate are clearly smaller in comparison with those for the 
cylinder. This is to be ascribed to the fact that in the case of the 
plate only the boundary-layer terms of the order e2 <<i (rectilinear 
flow) were neglected, whereas for the case of the cylinder there were, 
in addition, neglected terms of the order el <<1 that were related 
to the curvature. 


III. HYDRODYNAMIC AND THERMAL COMPARISON BETWEEN VERTICAL 
PLATE AND HORIZONTAL CYLINDER FOR FREE CONVECTION 
AND FOR PLATE IN PARALLEL FLOW 
1. Abstract 

Now that the theoretical boundary -layer solutions for the vertical 
plate and the horizontal cylinder for free convection have been obtained, 
it is of advantage to compare these two standard bodies of two- 
dimensional free flow with regard to the shape of the stream and the 
heat transfer. Similarly, a comparison of the boundary-layer develop- 
ment between the vertical plate for free convection and for the plate 
in parallel flow (Blasius solution of the boundary-layer equations) is 
of interest. The following features are characteristic for the three 
types of flow. 


^^The fact that the values of the heat-transfer coefficient ob- 
tained from the field measurements of Jodlbauer through integration of 
the temperature gradients at the wall over the cylinder perimeter lie, 
on the average, 4 percent higher than the theoretical values (see 
fig. 10) is due to the deviation between the theoretical and the actual 
conditions in the neighborhood of the upper stagnation point. 
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The houndary- layer thickness of the plate in parallel flov in- 
creases with the square root of the distance from the incidence 6dgej 
that of the vertical plate in free convection increases, however, as 
the fourth root of the height above the lower edge (cf. eqs. (69) and 
(71)). Both start with the thickness zero, which for the plate in 
parallel flow results in infinitely large velocity gradients, and for 
the free plate results in an infinitely large temperature drop and an 
infinitely large local heat -transfer coefficient at the lower edge. 

In the case of the cylinder, on the contrary, the boundary layer starts 
at the lower stagnation point of the horizontal cylinder with a finite 
thickness and thereafter with finite velocity and temperature gradients 
and a finite local heat-transfer coefficient (fig. 17). It increases 
according to a complex law (l/g(x), see eqs. (70) and (82)) and at the 
upper stagnation point reaches a theoretically Infinite thickness, with 
vanishing velocity and temperature drops normal to the surface (upward 
stream of warm air) and with a vanishing heat-transfer coefficient. 

In spite of these fundamental differences it is possible to set 
up a number of relations between the three cases. These show on the 
one hand a close hydrodynamic and thermal kinship which is due to the 
boundary- layer character of the differential equations underlying the 
theoretical solutions. On the other hand, they give practical view- 
points with regard to the application of rectangular plates or hori- 
zontal cylinders for the heat transfer. For this purpose it is first 
of all necessary to represent the Reynolds numbers of the characteristic 
length, of the boundary-layer thickness, and of the nondimens lonal local 
and average heat-transfer coefficients as functions of the Grashof 
numbers of the plate and cylinder, which in turn requires the introduc- 
tion of a boundary-layer thickness (flow discharge thickness) for the 
velocity profile of the free convection. A knowledge of the dependence 
of Re of the boundary-layer thickness on Gr of the plate, as well 
as of the azimuth and Gr of the cylinder, is in addition required 
for the later evaluation of the tests as regards the occurrence of 
turbulence (Part IV) for free^ convection at the plate and the cylinder. 


2. Bo\mdary-Layer Thickness for Free Convection 

It is first necessary to determine the magnitude to be associated 
with the boundary-layer thickness. For the velocity profile u(y) of 
the plate in parallel flow (U = maximum velocity), Prandtl and his 
coworkers, as is known, introduced the displacement thickness 5 , 
defined as 

5*- U = (U = u) dy (66) 

which is a measure of the decrease in the flow discharge as a result 
of the friction. 
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For the velocity profile of the free convection u(y), a flow 
discharge thickness 5 is introduced which is hydrodynamically equiv- 
alent to the displacement thickness, being like the latter character- 
istic of the development of the friction layer, which in this case 
takes up the entire flow; 

'^00 

5 • U = u dy (67) 

J 0 


or expressed in words: Through the flow discharge thickness 6 would 

flow the same fluid mass with the maximum velocity as actually flows 
with the total stream. 


In the nondimens ional velocity profile p'(q) (q = nondimens ional 
distance from the wall, p = stream function) which occurs in the theo- 
retical solution, the discharge flow thickness is denoted by qg, the 
maximum velocity by p'_^. From equation (67) there is then obtained: 

iTlcxX 


q.cP' 

^6 max 


p' (q)d<i = p(*) - p(o) 


and with the numerical values: 


P 


T 

max 


0.275 


p(oo) = 0.60^^ 


p(0) = 0 

the value of is obtained as 

qg = 2.18 (68) 

For orientation purposes^i, it is assumed that the velocity maximum 
lies at q = 0.95 and the point of inflection of the profile at 
q = 1.85. 


With the values of q^ for the plate and cylinder from equa- 
tions (72) and (82), respectively, there is obtained for the vertical 
plate, where h is the distance from the lower edge: 


^Extrapolation from last computed value p(6.0) = 0.5928. 
^^See continuous curves in figures 13 and 15. 
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I = 2.18 (69) 

and for the horizontal cylinder of diameter d at azimuth x 

I g(x)2V^ = 2.18 (70) 


3. Nondimens ional Representation of Theory for Vertical Plate 

(a) Reynolds numbers . - From the theoretical solution for the 
vertical plate given by Schmidt and Pohlhausen (ref. 21) the following 
expression is obtained according to their equations (23) , (24)^ and 
(31)^^^ since c^ = Grj^/4h^^ for the velocity u and the distance 

from the wall y: 

U = P’ (q) Gr^VS ^ = q ^2 ( 71 ) 

For the flow discharge thickness y = 6 and the maximum velocity 
u = U, there is then obtained 

U - X ^iax |=<!5V2 0rh-^A (72) 

or 

^ = 2V2 Piax 96 ^ - 2 Piiax (’5) 

and using the abbreviations = U*h/v and Rg = U*6/v and substi- 
tuting the numerical values for and q 

Rg = 1.695 GrjjV^ ( 74 ) 

Rj^ = 0.550 Grj^ ( 75 ) 

Elimination of Grj^ from the last two equations gives the following 

expression as the relation between Re of the flow thickness and the 
height of the plate for a vertical plate in free flow: 


^^The distance from the lower edge is denoted by h instead of x, 
the nondimens ional distance from the wall ^ by q^ the nondimens ional 
stream function ^ by p. 
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R 


5 


2.29 



(76) 


(Td) Nusselt characteristics. - For the local heat-transfer coef- 
ficient a(h) at the height h ^ahove the lower edge, the following 
equation holds (ref. 2l)*: 


a(h) 



(77) 


(t ’ = -0.508 is the value of the derivative of the nondimens ional 

temperature function t(q) for q = O) or, in nondimens ional form, the 
numerical values may he substituted: 

= 0.359 (78) 

The variation of the heat-transfer coefficient along the variable _ 
height h of a plate of total height H is most conveniently obtained 
by introducing a local heat-transfer coefficient which has been made 
nondimens ional by dividing through by H 

= Nujj(h) = 0.359 GrjjlA (79) 

The heat-transfer coefficient (fig. 18) at the lower edge _ of the plate 
with zero boundary-layer thickness is theoretically infinitely 
and then continuously drops with the reciprocal of the fourth root of 
the height. 

The mean heat-transfer coefficient corresponding to a plate of 
height H 

om = I a(h)dh (j) 

Jo 


is obtained from equation (77) as 


or, 

as 


nondimensionalized and with numerical values substituted (fig- 18), 

^ = Nujj = 0.479 GrjjV^ 


* Ext r apo lat ion . 


(81) 


MCA TM 1366 


39 


4. Nondimensional Representation of Theory for Horizontal Cylinder 

(a) Reynolds numbers . - The theoretical solution given previously 
for the tangential velocity u (eq. (55a)) and the distance n from 
the wall (eq. (55d)) for the azimuth x of a horizontal cylinder, 
with the values n = 6 and the maximum velocity u = U, may be trans- 
formed to 

U = ^ f(x)g(x)pi^ Gr//2 I = ^ Gr^-lA (82) 

For the upper stagnation point x = Jt, because g(rt) = 0 and 
f(jt) = 5.84, the boundary -layer thickness becomes infinite and the 
tangential velocity zero. The quantity flowing through U5 and there- 
fore the Re of the boundary layer U&/v, however, remain finite. To 
introduce Re for the characteristic length results in no simplication 
for the cylinder because, in addition to the characteristic length, the 
cylinder radius occurs as an additional characteristic length of the 
system. From equation (82), 

U5/v = qg Gr^l/4 (83) 

Substituting the numerical values and using Gr^ result in the follow- 
ing expression for Re of the boundary layer at azimuth x: 

Rq(x) = 0.357 f(x) Gr^V^ (94) 


and for Re of the boundary layer at the upper stagnation point 
(with f(«) = 5.84): 

Rg = 2.08 Gr^V^ (85) 

where Rg is written for R^(at). Although the theoretical solution 

for X = It is no longer valid^ equation (85) because of the small 
change of f(x) gives the correct relations in the region of the upper 
stagnation point (up to about 165^)^ so that for the sake of simplic- 
ity the formula will be used in the following discussion. 

(b) Nusselt characteristics . - Figure 18 shows the variation of 
the nondimensional local heat-transfer coefficient according to equa- 
tion (60) over the developed semicircumference of the cylinder. In 
the region of the lower stagnation point, it has a finite, practically 
constant value; and at the upper stagnation point, it drops rapidly to 
zero in accordance with the upstreaming warm air at that point. In the 
same figure the mean heat-transfer coefficient Nu^ for a cylinder 

is plotted according to equation (65). 
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5. Hydrodynamic Comparison 


(a) Between vertical pla~be and horizontal cylinder in free con - 
vection. - Although the boundary- layer growth for a vertical^plate and 
a horizontal cylinder, as already remarked (Part III, sec. l) and as 
seen from equations (72) and (82), differ fundamentally in character, 
it is nevertheless of interest to compare the flow condition (always 
characterized in what follows by Rg) at the upper edge of the plate 

and in the region of the upper stagnation point of the cylinder. The 
following three questions will here be considered: 

1. The behavior of the flow at the upper edge of the plate and at 
the upper stagnation point of the cylinder for the case that the height 
of the plate is equal to the cylinder diameter. 

2. The behavior of the flow at the upper edge of the plate and at 
the upper stagnation point of the cylinder for the case that the height 
of the plate is equal to the developed semicircumference of the cyl- 
inder, that is, the characteristic lengths are equal. 

3. The relation of the height of a plate whose flow condition at 
the upper edge is equal to that of the upper stagnation point of a 
cylinder to the diameter of the cylinder. 

Before each of these questions is considered, the ratio of the 
boundary-layer Re raunbers at the upper stagnation point of the cyl- 
inder Rg(z) and at the upper edge of the plate Rg(P) are written, 

with the aid of equations (74) and (85): 


where R 5 (z) and Gr^ refer to the cylinder, and Rg(P) and Gr^ 

refer to the plate. The coefficient 1.226 has the following signifi- 
cance (cf. eqs. (85) and (83) for the cylinder, eqs. (74) and (73) for 

the plate) : 


R5(Z)/Rs(P) = 1.226 


( 86 ) 



(87) 


1. For equal Gr of plate and cylinder,' Gr^^^ = Grj^, equation (86) 


gives 


Rs(z) = 1 . 226 -R 5 (P) 


( 88 ) 
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that is, for equal Gr of plate and cylinder, for example, for a plate 
of the height of the. cylinder diameter h = d under otherwise equal 
conditions (with regard to temperatures and materials). Re at the 
upper stagnation point of the cylinder is 22 percent greater than that 
of the upper edge of the plate (fig. 19). This result states nothing 
about the hydrodynamic relation between the two, since the comparison 
refers to different characteristic lengths. It is useful only for the 
rapid comparison of plate and cylinder for the usual values of Gr. 

It is then known that for *^qual Gr the cylinder possesses the more 
developed flow, so that the flow at the cylinder may be turbulent while 
that of the plate is still laminar. 

2. For equal characteristic length h = jt/2d there is obtained 
from equation (86) 

R5(Z)/Rg(P) = 1.226 (2/n)V^ = 0.874 (89) 

The boundary-layer Re at the upper stagnation point of a cylinder is 
thus 13 percent smaller than that of a plate of the height of the 
developed semicircumference of the cylinder (fig. 19). This retar- 
dation of the boundary-layfr development as compared with the vertical 
plate is due to the fact that the boundary layer at the lower stagnation 
point of the cylinder already has a finite thickness, whereas in the 
case of the plate the thickness must increase from zero. 

3. For equal boundary-layer Re values for the plate and cylinder, 
equation (86) gives 

GrgjGr^ = 1.226"^ = 0.441 (90) 

At the upper stagnation point of the cylinder, there is therefore the 
same flow condition as at the upper edge of the plate if Gr of the 
cylinder is 0.44 times that of the plate. This permits a rapid con- 
version of the flow data for the plate into those for the cylinder and 
vice versa. If, for example, it is known from heat-transfer experiments 
on a plate that the departure from the laminar GrV^ law occurs at 
(Grh)kr = 8 ^ 108 , then it can be concluded without experiment from equa- 
tion (90) that this must be the case for the horizontal cylinder for 
(Grd)kr = 3 . 5 x 10 ®. If the plate and cylinder are under otherwise equal 

conditions, it is possible from equation (90) to derive a simple rela- 
tion for the height of a plate hg^ which is hydrodynamically equivalent 

to a cylinder: 

d.3;hg3 = 0.441 hg:d = 1.311 


(91) 
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In the neighborhood of the upper stagnation point of a horizontal cyl- 
inder, the same flow condition prevails (Rg) as under otherwise equal 

conditions at a vertical plate of height 1.31 times that of the cylin- 
der diameter (fig. 19). As was to be concluded from the answers to 
questions 1 and 2, this height must lie between the cylinder diameter 
and the developed semicircumference. 

The answer to the three questions on the flow condition at a 
cylinder and at plates of various heights is shown in figure 19 (Rg 

of the cylinder set equal to l) . 

(b) Hydrodynamic comparison between a vertical plate in free con - 
vection and a plate in a parallel flow . - The relation (76) between Re 
of the boundary-layer thickness and that of the height of the plate 
(characteristic length) vertical plate in free convection 

is strongly analogous to the corresponding relation (92) between Re 
of the boundary-layer thickness Rg* (referred to the displacement 

thickness 6*) and the characteristic length R^^ for a plate in a 
stream parallel to its plane (Blasius solution of the boundary-layer 
equations) : 

Rg = 2.29 • R^^/^ (76) II Rg* = 1.73 R^^/^ (92) 

In spite of the already mentioned entirely different rate of growth of 
the boundary-layer thickness, in the case of the parallel flow as the 
square root of the distance from the incident edge and in the case of 
the free convection as the fourth root of distance from the lower edge, 
the development of the Reynolds number of the boundary-layer thickness 
as a function of the Reynolds number of the characteristic length in 
the case of these two very different types of flow follows the same 
exponential law and with approximately equal magnitude. 


6. Thermal Comparison between Vertical Plate and Horizontal Cylinder 

(a) Mean heat-transfer coefficients . - For equal Gr of plate and 
cylinder (Gr^ = Grg^), there is obtained from equations (8l) and (65) the 

following ratio of the nondimens ional heat-transfer coefficients of the 
horizontal cylinder Nug^ and of the vertical plate Nug 

Nu^/Nug = 0.777 (93) 

Under otherwise equal conditions, therefore, the mean heat-transfer 
coefficient of a cylinder is only 78 percent of that of a plate of the 
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height of the cylinder diameter (fig. 18). In agreement with this, the 
comparison given hy M. Jakob and W. Linke (ref. 19) of convection tests 
on vertical plates and horizontal cylinders in a single Nu-Gr diagram 
(with H and d denoting the length dimension) gives a lower position 
of the cylinder tests. The modified determination of the critical num- 
ber, which follows from these results, will be discussed later (Part IV, 
sec. 5). 

For equal mean heat-transfer coefficient of plate and cylinder 

under otherwise equal conditions, there is obtained from equations (8l) 
and (65) 

0.479 = 0.372 d“V^ H/d = 2.76 (94) 

The mean heat-transfer coefficient of a cylinder is therefore, under 
otherwise equal conditions, equal to that of a plate of height 2.76 
times that of the cylinder diameter. On account of the small decrease 
of the local heat-transfer coefficient at large plate heights, this 
factor of 2.76 is not very sharply determined. Thus Jodlbauer (ref. 29) 
found its possible, without considering the previously given theoretical 
solution for the cylinder, to represent the mean heat-transfer coef- 
ficients for cylinders from tests by Koch by the formulas for the ver- 
tical plate if H is replaced by 2d (instead of theoretically by 
2.76d). As a matter of fact, the values thus determined for the 

cylinder lie only 8.3 percent above the values given by the theory 
according to equation (65), as can be seen by replacing H by 2d in 
equation (81) . 

(b) Total heat transfer . - The total heat transfer Q(z) of a 
cylinder of diameter d along the entire circumference and the total 
heat transfer Q(P) along the twc sides of a plate of the height H, 
under otherwise equal conditions, are in the ratio 

Q(Z) _ 

Q(P) ■ Om(P)2H 2 Nuh(P) 

and from eqimtions (8l) and (65), 

Q(Z)/Q(P) = 1.221 Gr^l/VGr^V^: (95) 


where the significance of the factor 1.221 is as 
(80), (65), and (59)): 


1.221 


jt 0 . 37 2 _ at 
2 0.479 “ 2 


2V^g(-to') 


4 ^ 
3 2 


(-to’) 


follows (cf. eqs. (8l), 


3otg 


(96) 
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A comparison of this factor 1.221 with the factor 1.226 occurring in 
the ratio of the houndary-layer Re of cylinder and plate (eq. (86)) 
according to equation (87) shows that hoth are identical, since accord- 
ing to the theoretical solution, the following relation holds as an 
immediate consequence of the validity of the heat-transfer equation: 

3jtg = f(jt) (97) 

The small difference of l/2 percent is explained hy the fact that (as 
mentioned in connection with eq. (64)) g was computed not according 
to equation (97) hut was obtained hy planimetering. 

For equal Gr^ and Grjj values in equations (86) and (95)^ the 
following is thus obtained: 

Q(Z)/Q(P) = R5(Z)/R5(P) (98) 

that is^ the total heat transfer of a horizontal cylinder and that of 
the two sides of a vertical plate under otherwise equal conditions are 
in the same ratio as the boundary -layer Re values at the upper stag- 
nation point of the cylinder and at the upper edge of the plate^ 
respectively. 

On the basis of relation (98) it is now possible, without further 
computation, to answer the three questions on the total heat transfer 
of a cylinder and plates of different height, analogous to the questions 
referring to the hydrodynamic comparison (see fig. 19; the heat transfer 
of the cylinder has been set equal to l). 

Corresponding to equation (88), 

Q(Z) = 1.22 Q(P) (99) 

The total heat transfer of a cylinder is thus 22 percent greater than 
that for the two sides of a plate of a height equal to the cylinder 
diameter under otherwise equal conditions (fig. 18). 

Corresponding to equation (89), 

Q(Z) = 0.87 Q(P) (100) 

The total heat transfer of a cylinder is 13 percent smaller than that 
of the two sides of a plate of height equal to the developed semi- 
circumference of the cylinder (equal characteristic length) under other- 
wise equal conditions. 


NACA TM 1366 


45 


Corresponding to equation (91 ) for 

Q(Z) = Q(P) 

there is obtained the height of a plate H- thermally equivalent to 
a cylinder: 


% : d = 1.31 (101) 

The total heat transfer of a cylinder under otherwise equal conditions 
is therefore as great as for the two sides of a vertical plate of 
height 1.31 times that of the cylinder diameter. In summary it may 
therefore he stated that a horizontal cylinder and a vertical plate^ 
which transfers heat on both its sides, of a height equal 1.31 times 
the cylinder diameter are equivalent both thermally (with reference to 
the total heat transfer) and hydro dynamic ally (with reference to the 
boundary-layer development at the upper stagnation point and at the 
upper edge, respectively). 


IV. THE OCCUKEIENCE OF TURBULENCE IN FEEE FLOW ABOUT A HORIZONTAL 
CYLINDER AND ALONG A VERTICAL PLATE 
1. Introductory Observations 

For the determination of the transition from laminar to turbulent 
flow in the case of free convection, tests were carried out on a verti- 
cal plate and a horizontal cylinder of sufficient size. Two main 
reasons underlay the investigation. In the first place, it was desired 
to determine the upper limit of the validity of the laminar boundary- 
layer theories for plate and cylinder and therefore the upper limit of 
applicability of the heat-transfer formulas developed from these theo- 
ries for practical application. The second reason was of a more 
theoretical nature. The numerous turbulence investigations, both ex- 
perimental and theoretical, undertaken in recent decades have been 
concerned almost exclusively either with the flow between two parallel 
walls with linear velocity distribution (Couette flow) or with the flow 
in pipes, channels, about rotating cylinders, along plates, cylinders, 
or other resistance bodies which have in common a velocity profile which 
rises uniformly from the value zero at the wall up to a maximum value. 
Differing essentially from these profiles are evidently velocity pro- 
files for which the velocity rises from the value zero at the wall to 
a maximum value and then, after passing through a point of inflection, 
again drops to the value zero at a large distance from the wall. Such 
profiles are typical in the case of heat transfer by free convection. 
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Some further observations on the setting up of turbulence in free con- 
vection will be made later in connection with the discussion of results 

(sec. 5). 


2. Test Procedure 


For the investigation, the schlieren method developed "byE. Schmidt 
was applied, which as an optical method provides the possibility of an 
instantaneous view of the entire field, and is therefore of advantage 
for determining the critical number. 


The schlieren method has been described in detail by E. Schmidt^ 
(ref. 22) so that only the essential points will be given here. It is 
based on the deflection suffered by a light ray in passing through a 
field with density stratification (i.e., variation of the index of 
refraction) normal to its direction of propagation toward the colder 
layer. In the test, a parallel light beam is allowed to pass tangen- 
tially along the surface of the heated body. The rays in the neighbor- 
hood of the wall, because of the maximum temperature ^op at that place, 
are most strongly deflected from their initial direction; those farther 
away from the wall are less deflected, while those rays which are out- 
side the temperature field undergo no deflection. On a screen which 
is set up at a sufficient distance behind the heated body, the rays 
near the wall appear the farthest from the wall. The following picture 
is obtained: the heated body throws a completely dark shadow which, as 
compared with that of the cold body, is increased by the thickness of 
the temperature boundary layer. This shadow is bounded by a fi^st 
"interior caustic curve" adjacent to which is formed a medium-bright 
region (rays which have passed through the temperature field), which in 
turn is bounded toward the outside by a second "external caustic curve 
originating from the tangential rays in the neighborhood of the heated 
surface. Still farther toward the outside is the uniformly bright field 
of the undisturbed illiunination. 


For an accurate evaluation, the method is suitable only for bodies 
having one of its length dimensions (two-dimensional problem), so that 
the effect of the temperature fields at the ends on the path of the 
light ray is small compared with the deflection undergone in passi^ 
through the distance along the length. On the other hand, the len^h 
should be chosen only large enough so that the light ray at the end of 
the "body still is approximately at the same distance from the wal as 
at the initial point of the body. 
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3. Test Setup and Procedure^® 

(a) Plate- - The plate was 100 by 100 centimeters and 1-0 centi- 
meter thick. It consisted of two zinc sheet plates 1 millimeter thick 
which by interposing a frame of 8-millimeter U-brass were soldered onto 
each other. In the interspace there was placed a heating coll, insu- 
lated by asbestos boards, which consisted of 36 meters of chrome-nickel 
wire of 1-millimeter diameter with a total resistance of 51 ohms. In 
air of room temperature, the plate could be kept at 100° C surface tem- 
perature with 5.6 amperes. The plate was suspended on two bicycle 
wheel spokes which were attached to the upper strip of the frame. 

Through the latter the current was also conducted to the Interior. For 
measuring the surface temperature a silver-constantan thermocouple was 
used which was flatly soldered on about 5 centimeters from the upper 
rim of the plate and the wires of which extended a few more centimeters 
quite close to the plate. This simple temperature measurement could be 
applied because extreme accuracy was not required. Figure 20 shows 
the plate suspended with the current -supplying wires, and the thermo- 
couple and two pendulimis each with two spheres which served to determine 
the vertical and made possible the mutual comparison of the different 
photographs . 

(b) Cylinder. - The cylinder consisted of two layers of 0.7- 
millimeter-thick zinc sheet which was rolled over two bicycle wheel 
rims and soldered to their upper sheathing strip. Between the sheets 
there was placed the heating coil insulated by asbestos board, the coil 
being of 54 meters of chrome-nickel wire of 1-millimeter diameter with 
a total resistance of 77 ohms, which brought the cylinder in air of 
room temperature to 100° C surface temperature with a ciirrent of 4.1 
amperes. The cylinder was suspended on two bicycle wheel spokes. The 
current leads and thermocouple were similarly located on the upper 
sheathing strip. The mean length of the cylinder was 100.2 centimeters, 
and its mean outer diameter 58.45 centimeters (fluctuating between 
58.35 in the horizontal and 58.75 over the soldered seam). The frontal 
surfaces were closed off with asbestos board. Figure 21 shows a view 
with the suspension, current leads, and the two pendulums. 

(c) Setup. - Since it is very difficult to produce a parallel 
bundle of rays of the diameter of the body here employed by means of 
lenses, such a beam was replaced by the light from a strongly screened 
(as a rule of 3-mm aperture) arc lamp 32 meters from the heated body. 

The distance between the body and the screen was 8 meters. The appa- 
ratus was set up on the lajrgest floor, of 42-meter length, of the 


^^Acknowledgment is made to Herr Dr. R. Weise and Dr. H. Kurzweg, 
Leipzig. For the test setup, acknowledgment is made to Master Mechanic 
G. Hentsch, Leipzig. 


48 


WACA TM 1366 


Institute. Because of the window in one of the walls and the connec- 
tion to the staircase halls, it was not possible to exclude entirely 
currents- of the room air so that only at high temperature differences 
(above 50° C) could useful schlieren photographs be obtained. The 
existing dissymmetry of the critical points on the two sides of the 
bodies, both for the plate and cylinder (see table V) was due to a 
unilateral convection flow from the window wall. The horizontal 
adjustment of the test bodies was effected by means of a water balance; 
the adjustment parallel to the light cone, by. measuring the shadow of 
the cold body. 

(d) Photography^^ . - The schlieren pictures were obtained by two 
methods. In one method the shadow pictures were taken on the screen 
by means of a lens and camera (13x18 cm, Zeiss-Tessar 1:2.7). In spite 
of highly sensitive plates (Agfa-Superpan 16/ 10 DIN), exposure times 
of 5 to 25 seconds were required on account of the low brightness, be- 
cause of the strong screening and large distance from the arc lamp. 

For determining a time mean value of the fluctuating critical points, 
these time pictures were, however, particularly well suited. The 
photograph was taken either somewhat from the side (plate) or, in order 
to avoid distortions, exactly central (partly in the case of the cyl- 
inder) ; the camera was set up in the shadow of the cylinder so that the 
legs of the tripod threw an additional shadow (fig. 23) . The second 
method consisted of a direct illiimination of one or several photoplates 
simultaneously (Agfa-Isochron 13/l0 DIN, 10x15 cm), which were brought 
to the position of the screen. The greater brightness obtained in this 
way made possible instantaneous photographs (through light flashes from 
the arc lamp), which in connection with their greater scale provided 
a view in instantaneous detail of the transition from laminar to turbu- 
lent flow. 


4. Evaluation of the Schlieren Pictures 

Figures 22 and 23 show time exposures of the plate and cylinder 
taken by means of lens and camera. In the center the magnified shadow 
of the body by the thickness of the temperature boundary layer, the 
"inner caustic curve" surrounding it, and further toward the outside 
the "outer caustic ciirve" arising from the rays near the wall. The 
dissymmetry of the lower edge of the plate is due to the disturbances 
of the room air that were mentioned previously. The upper edge of the 
plate is immediately below the frame of the picture in the "necking" 
of the central shadow. In the upper stagnation point of the cylinder 
are seen the shadows of the suspension wires, current leads, and thermo- 
couple wires; below are seen the three tripod legs for the central 


^"^For valuable advice in connection with the photography, ac- 
knowledgment is made to Herr Dr. A. Naumann, Leipzig. 
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photograph. Figure 24 shows for simultaneously obtained instantaneous 
photographs of the left^® side of the plate (in the case of direct 
illumination about 25 cm from the lower edge starting with small inter- 
spaces of about 1 cm), a region of about 50-centimeter height: at 24(a) 
and the lower part of 24(b), laminar; in the upper part of, 25(b), a 
transition; 24(c) and 24(d), turbulent. Figure 25 shows three simul- 
taneously taken instantaneous pictures of the cylinder with direct 
illumination: 25(a) at the lower stagnation point, laminar; 25(b) at 
the left azimuth of 120°, turbulent; 25(c) at the upper stagnation 
point, warm air stream. 

As shown by figures 22 to 25, the outer caustic c\irve of the rays 
near the wall is particularly suitable for determining the critical 
location. If this line is broken up, then turbulence certainly exists; 
whereas the flow regions farther from the wall are more easily disrupted 
through external disturbances. In the laminar (lower) region, the 
caustic curve is sharp and is completely stationary; in the turbulent 
(upper) region it is torn into several pieces and in disordered motion. 
In the time photographs, this means complete washing out. The critical 
position itself fluctuates upwards and downwards, at one time very 
rapidly and at another slowly. For evaluating the time photographs 
there was taken as the "time mean value of the critical position" the 
end of the region in which the existence of an outer caustic curve 
through differences in brightness could still be determined outwards 
and inwards (arrows). Above the critical position the heated layer 
rises as a whole outwards as broad, uniformly bright strips. Comparison 
of the photographic observation thus determined with the mean critical 
position estimated from naked eye observation over a period of time 
showed good agreement and therefore the justification of this procedure. 

The critical positions obtained by this procedure are collected 
together with the remaining test values in table V. Each critical point 
is the mean value of four determinations two of which were made with 
direct illumination, and two with a magnified projection picture. The 
Gr values are formed in the case of the cylinder with the diam'eter, in 
the case of the plate with the critical height as characteristic length. 
By means of equation (74) for the plate and, with f(xjj^p), equation (84) 

for the cylinder, the critical Reynolds numbers of the boundary appear- 
ing in the last column were computed. Their mean values of 303 for the 
plate and 285 for the cylinder are weighted according to the different 
illimination times of the individual pictures. From the mean value 

= 303 the value 1.0x10^ is obtained for the mean critical Gr 

for the plate. 


^®0n the right of the figure; the pictures with direct illimination 
are laterally interchanged. 
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5. Discussion of Results 

The mean value of 303 for the critical Reynolds number of the 
boundary later for the plate and 285 for the cylinder can be regarded 
in agreement, in view of the accuracy of determination of a critical 
number in general and the small number of the schlieren pictures and 
disturbances of the room air in particular. This is understandable 
from the fact that the velocity profiles of both flow processes at any 
distance from the lower edge or lower stagnation point are affine to 
each other. The development along the initial stretch occirrs in a 
different manner, but this should have only a higher-order effect. As 
was to be expected, the critical number of the free profile Rg “ 300 

lies considerably below the critical number of the profile in the case 
of forced flow (plate), which, referred to the displacement thickness 
of Burgers and Hansen, amounts to 950 in the wind tunnel, and accord- 
ing to more recent tests in the Gottingen water tunnel for very quiet 
approach, it amounts to about 1400. It would be of interest to know 
whether a theoretical stability investigation for the free profile, the 
analog of the familiar computations in the forced-flow case, would have 
to be conducted in order to give this increased instability. 

With the determination of the critical number for the vertical 
plate and the horizontal cylinder and reduction to the critical Reynolds 
number of the boundary layer, the data of other investigators will be 
compared with regard to the transition from laminar to turbulent flow. 

E. Schmidt (ref. 26) notes that according to tests in air, the flow 
at about 100° C plate temperature remains laminar up to a plate depth 
of about 50 centimeters, corresponding to a critical Gr of 8><10®, in 
good agreement with the values obtained herein. 

M. Jakob and W. Linke (ref. 19) determined from a combination of 
the measurements of various investigators on vertical plates, vertical 
and horizontal cylinders, block, and sphere, as the transition point 
from the GrV^ to the GrV^ law, the critical number of 
(Gr.Pr)i^j. 3.0 x 10'7, that is, for diatomic gases Gr^^-j, » 4.0x10". The 

difference as compared with the present determinations should in large 
part be -ascribed to the fact that hydrodynamically unequivalent bodies 
(Reynolds number of the boundary, see Part III, sec. 5a, eq. (88)) were 
here treated together. If account is taken only of the tests on the 
vertical plates and vertical cylinders (Part III, sec. 6a), which for 
large cylinder diameters and large Gr (thin boundary layer) are at 
least approximately equivalent hydrodynamically, there is obtained a 
critical Gr of 4x10®, which is therefore quite close to the present 
determinations . 
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E. Griffith and A. H. Davis (ref. 17) investigated the local 
distrihution of the heat transfer on a 270-centimeter -high wall, which 
consisted of 25 individual elements. They obtained a minimum value of 
the heat transfer at about 38-centimeter height (their own determina- 
tion) which they regarded as the transition from the laminar to the 
turbulent condition. This agrees in order of magnitude with the pres- 
ent values (columns between the individual elements), where it shoiiLd be 
further clarified which point of their curve is to be regarded as the 
"critical point" proper. 

6. Travel of the Critical Azimuth for the Horizontal Cylinder 

On the basis of the given theoretical solution of the flow about 
the horizontal cylinder and the point of transition experimentally 
determined for a Gr “ 109, the travel of the critical azimuth for a 
variable Gr can be computed if it is assumed that the turbulence 
always occurs if the boundary-layer Re has the same value of 285. 

This reasonable assumption will correspond all the more to fact since, 
as has just been shown, even for the vertical plate with a quite differ- 
ent boundary- layer development, the turbulence is set up at the same 
critical number of the boundary layer. 

From equation (84), the following relation is obtained, with the 

value Rc, = 285: 

°kr 


which is evaluated in table VI and plotted on semi logarithmic grid in 
figure 26. In accordance with the basic considerations, for a Gr = 10 
the start of turbulence occurs for an azimuth of about 120°. Complete- 
ly laminar flow about the cylinder, for extension of the theory up to 
the upper stagnation point, occurs for Gr<3><108. The transition poin 
reaches the equator for Gr = 3X10^ . 

For greater clarity, the Gr will be expressed for cylinders of 
varying diameters, which are in air of 20° C with a surface temperature 
of 100° C. From the definition of Gr, for 0 = 80 (°C), 

6 = 1/293 (°C"^), v(i00°) = 0.231 (cm^/sec), the following relation 
obtains : 


which with the aid of the critical Gr^ values from equation (l02) has 

similarly been evaluated in table VI and plotted in figure 27 . Corre- 
sponding to the basic considerations, the critical azimuth for the 



( 102 ) 



(103) 
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relations considered for a cylinder of 60-centimeter diameter is 120°. 
For cylinders "below 41 centimeters, no turbulence occurs at all. The 
transition point reaches the equator for a cylinder of 84-centimeter 
diameter. 

For practical applications, the result is thus obtained that for 
heated horizontal conducting pipes up to 40 centimeters in diameter and 
surface temperatures of 100° C in room-temperature air the previously 
given formulas are valid for the heat transfer in air (in particular, 
eqs. (60) and (65)). For accurate computation, the nondimens ional 
formulation of figure 26 is naturally to be applied. The formulas can, 
in the absence of other suitable data, still be applied up diameters of 
60 centimeters, because two thirds of the cylinder periphery is still 
in laminar flow. There should not likewise be much of a change for 
considerably higher surface temperatures (200° or 300° C). 

The experimental determination of the variation of the critical 
azimuth for the cylinder with the Or, through varying the temperature, 
remained unsuccessful. For temperature difference up to about 50° C 
the outer caustic curve is still not sufficiently far from the central 
shadow in order to be perceived separately from the latter, as is re- 
quired for determining the critical munber. For this, the distance 
between the cylinder and the screen would have to be considerably in- 
creased, which is not possible on account of space requirements. Going 
beyond the usually employed temperature differences of 80° C, because 
of the strong increase in v, no longer gives an increase in Gr and 
moreover gives increasingly stronger deviations from the assumption of 
moderate temperature differences, which is at the basis of the theory. 
The travel of the critical azimuth with Gr is best determined by 
varying the cylinder diameter, as follows from figure 27. 


V. SUMMARY 

From the numerous tests already available on the heat transfer 
from horizontal pipes and wires in diatomic gases, the dependence of 
the nondimensional heat -transfer number Nu on the Grashof number Gr 
and the temperature coefficient Te, which enters as a further non- 
dimensional factor at large temperature differences, is determined. The 
effect of Te on the heat transfer is quantitatively determined for the 
first time. It is particularly large in the region of small Gr (lO"'^ 
to lO), where on the average a decrease in Nu by 22 percent for in- 
crease of Te from Te = 0 to Te = 1 is obtained. The experimentally 
found dependence of Nu on Gr, Pr, and Te can then obtain a qualita- 
tive theoretical explanation. 
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For the region of Gr (about 10^ to 3^10®) in which the heat trans- 
fer is limited to a thin (in comparison with the cylinder diameter)^ 
heated layer with laminar flow^ the velocity and temperature fields and 
the heat transfer are quantitatively computed from the boundary-layer 
differential equations without any additional empirical values; these 
computations are found in good agreement with the available measure- 
ments. In particular^ the one-fourth power law of the heat transfer is 
obtained theoretically as Nu = 0.37 GrV^. 


The flow and heat -transfer relations thus computed (Re of the 
boundary-layer flow, variation of the local heat -transfer coefficient 
along the cylinder periphery, mean heat -transfer coefficient, and total 
heat given off by the cylinder) are compared with the already known 
relations for the vertical plate. Among other results, the depth of a 
vertical plate is detemined which for free convection shows equal flow 
condition and equal total heat transfer at a given cylinder, so that a 
simple computation is made possible for converting the heat-transfer 
data for a plate to those of a cylinder and conversely. 


In order to know the upper limit of validity of the laminar -flow 
and heat-transfer computations and the laminar heat-transfer formulas, 
the start of turbulence was determined by schlieren photographs on a 
vertical plate and a horizontal cylinder of sufficient size. This 
occurs for different values of Gr, namely, Gr = 1.0^10® for the plate 
and Gr = 3.5^10® for the cylinder, but for equal Re of the boundary- 
layer flow, namely, for Re * 300. In this way, the critical Reynolds 
number of the velocity profile of the free flow was determined for the 
first time; as compared with other velocity profiles, this profile is 
marked by the presence of a maximum value and a point of inflection. 
Likewise computed with the Grashof number was the * critical azimuth* on 
the cylinder at which the laminar flow passes into the turbulent flow. 
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TABLE I. - SUMMARY OF TESTS BY DIFFEREMT INVESTIGATCRS 


Investigator 

Ayrton, 

Kilgour 

Langmuir 

Bijlevelt 

Kennelly, 

Wright, 

Bylevelt 

Petavel 

Wamsler 

Koch 

Substance 

9 Pt wires 

5 Pt wires 

6 wires: 

Ta, Pt, Fe, 
Cu, Ag, Ni 

3 Cu wires 

1 Pt wire 

Pipes: 

6 wrought 
Fe; 6 Cu; 
1 cast Fe 

4 steel 
pipes 

Cylinder diam . , 
cm 

0.0031- 

0.0356 

0.00404- 

0.0510 

0.0043- 

0.1000 

0.01143- 

0.06907 

0.1106 

2.05-8.9 

5.9 

5.9 

1.4- 

10.05 

Temperature of 
body, t^, °C 

40-300 

227-1027 

46-239 



--- 

27.6-188.6 



Excess tempera- 
ture, 0, °C 

— 

— 


15-180 

200-1000 

36-243 

13-174.4 

Room tempera- 
ture, t„, °C 

10.5-15.9 

27 

17-22 

20 (approx) 

16 

10-29 

14.2-22 

Medium 

Air 

Air 

Air 

Air 

Air, Hg, 
°2 

Air 

Air 

Pressure, 
mm Hg 

750 

750 

750 

120-1900 

0.12- 
160 atm 

715 

715.7-722.3 

Cylinder length, 
cm 

32.5 

100 

35 

150 (approx) 

45 

300 

138-198 

Surrounding 

space 

Horizontal 
pipe; 
length, 
32.5 cm; 
diam. 

5.08 cm 

Base, lOOx 
15 cm2; 
height, 

30 cm 

Height, 
300 cm 

Vertical 
tank: 
length, 
152 cm; 
height, 
660 cm 

Horizontal 

pipe: 

length, 

45 cm; 
diam. 

2.06 cm 

Base, 7.35 
sq meters; 
height, 

210 cm 

Height, 
400 cm 

Ratio of space 
height to 
cylinder diam. 

1640: 143 

7430: 588 

70,000:3000 

58,200:9640 

18.6 

102:24 

286:40 

Te number 

0.30-1.00 

0.67-3.33 

0.097-0.758 

0.052-0.61 

0.69-3.46 

0.13-0.82 

0.045- 

0.606 

Relative 
decrease of 
Nu for 
ATe = 1, 
percent 

7.6-16.8 

7. 2-2. 2 

8.4-45.6 

25.2 

Irregular 

None 

None 

Mean value 
of Nu for 
ATe = 1, 
percent 

13.2 

Systematic 

25.6 

25.2 





Figure 

1 

2 

3 

4 

5 

6 

7 


CH-8 


TABLE II. - VALUES FOR THE MOST PROBABLE HEAT -TRANSFER LAW IN DIATOMIC GASES 


[Nu as function of Gr with Te as parameter for Pr = 0.74. ] 



log Gr 
Gr 

-4 

10"^ 

-3 

10"^ 

-2 

10"^ 

-1 

10"^ 

0 

1 

1 

loi 

2 

10^ 

3 

10^ 

4 

10“^ 

5 

10^ 

CD 

CD o 


GO 

00 o 

Te=0 

log Nu 
Nu 

0.685-1 

.484 

0.716-1 

.520 

0.787-1 

.612 

0.908-1 

.809 

0.040 

1.10 

0.176 

1.50 

0.321 

2.18 

0.476 

2.99 

0.650 

4.47 

0.850 

7.08 

1.094 

12.4 

1.344 

22.1 

1.594 

39.3 

Te= 

0.65 

log Nu 
Nu 

0.618-1 

.415 

0.649-1 

.446 

0.719-1 

.524 

0.838-1 

.688 

0.985-1 

.966 

0.134 

1.86 

0.290 

1.95 

0.458 

2.87 

0.645 

4.42 

Same as for Te = 0 


TABLE III. - NORMED AZIMUTH FUNCTIONS F(x) AND g'(x) AS SOLUTION OF DIFFE^NTIAL 
EQUATION SYSTEM (28) WITH BOUNDARY CONDITIONS ACCORDING TO EQUATION (29) 


[f^ solution of eq. (33l); F^^ solutions of eq. (33Il).] 


X 

0° 

30° 

60° 

90° 

120° 

150° 

165° 

180° 

F' 

F” 

G 

0.0007 

(-0.23) 

1.000 

0.518 

.456 

.989 

1.025 

1.018 

.945 

1.504 

1.504 

.873 

1.936 

1.941 

.765 

2.298 

2.352 

.602 

2.433 

(2.55) 

.473 

2.560 

2.763 

0 


TABLE IV. - FINAL AZIMUTH FUNCTIONS f(x) AND g(x) AS SOLUTION OF DIFFERENTIAL 
EQUATION SYSTEM (24) WITH BOUNDARY CONDITIONS OF EQUATIONS (26a) TO (26c) 


X 

0” 

30° 

60° 

90° 

120° 

150° 

165° 

180° 

f 

0 

1.181 

2.337 

3.430 

4.41 

5.24 

5.55 

5.84 

g 

0.760 

.752 

.718 

.664 

.581 

.458 

.360 

0 

f’-g 

0 

.887 

1.678 

2.274 

2.565 

2.397 

1.995 

0 
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TABLE V. - EVALUATION OF SCHLIEREN PHOTOGRAPHS FOR OCCURRENCE 

OF TURBULENCE 


(a) Vertical plate. 


Number 

b, 

mm 

°c 

Exposure 

time^ 

sec 

°c 

Critical height^ ^kr^ 

hkr 

^5kr 

Left 

Right 

Middle 

24 

748 

16.5 

20 

78 

76.0 

55.9 

66.0 

13.42 xlO® 

324 

25 

748 

16.5 

20 

113 

53.8 

47.3 

50.6 

6.83 xlO^ 

274 

26 

748 

16.5 

25 

100 

66.2 

60.8 

63.5 

13.13 xlO® 

323 

27 

748 

16.5 

20 

100 

60.9 

47.3 

54.1 

8.12 xlO® 

286 


Mean : 303 


(b) Horizontal cylinder. Diameter, d, 58.45 centimeters. 


Number 

b, 

boo, 

Exposure 

time^ 


Critical azimuth, x^^ 

G’^hkr 


^6kr 


mm 

°C 

Ori 

Left 

Right 

Middle 







sec 








6 

760 

19 

10 

102 

125.5° 

o 

123.5 

124.5° 

10.26 X 10® 

4.55 

291 

10 

760 

19 

20 

102 

125.8° 

110.0° 

117.9° 

10.26 X 10® 

4.35 

278 

11 

760 

19 

5 

102 

120.9° 

116.0° 

118.5° 

10.26 X 10® 

4.37 

279 

12 

760 

19 

20 

102 

126.7° 

117.6° 

122.1° 

10.26 xlO® 

4.48 

286 

17 

752 

18 

20 

102 

123.8° 

115.0° 

119.4° 

10.23 X 10® 

4.39 

281 

18 

752 

18 

5 

102 

144.5° 

129.5° 

137.0° 

10.23 xlO® 

4.90 

313 


Mean; 285 


TABLE VI. - TRAVEL OF TURBULENCE TRANSITION POINT (CRITICAL AZIMUTH 
WITH GRASHOF NUMBER OF CYLINDER Gr^ OR WITH CYLINDER DIAMETER d 
[t^, 100° C; t„, 20° C; b, 760 mm.] 


^kr 

180° 

165° 

150° 

120° 

o 

o 

60° 

30° 

15° 

o 

o 

H 

5° 

Gr j 

3.52 

4.30 

5.42 

1.07 

2.96 

1.37 

2.10 

3.35 

1.69 

2.72 

d 

xlo8 

xlo8 

xioQ 

XlO® 

XlO® 

X lObO 

xiobb 

xl0^8 

xl0b3 

xlOl^ 

d, cm 

41.3 

44.1 

47.7 

59.9 

84.0 

140 

347 

890 

1500 

3790 
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Figure 1. - Convection tests of Ayrton and Kilgour on nine platinum wires. 
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Figure 2. - Convection tests of Langmuir on five platinum wires. 
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Figure 3. - Convection tests of Bijlevelt on six wires of tantalum, platinum 


iron, copper, silver, and nickel. 



Figure 4. - Convection tests of Kennelly, Wright, and Bijlevelt on three 
copper wires. 
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Figure 5. - Convection tests of Petavel on one platinum wire. 



Figure 6. - Convection tests of Wamsler on one copper, one cast iron, and six 
wrought iron pipes. 
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Pr-O.T<* 

1 Ayrton and Kilgour 

2 Bijlevelt 

3 Kennelly et ai. 
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Figure 8. - Comparison of test results of different investigators for Te = 0 
(extrapolated). Limiting law of small temperature differences. 
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Figure 9. - Comparison of test results of various investigators 
for Te = 0.65. 
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Figure 10. - Summary of all test results with diatomic gases (Pr = 0.74) with 
the adjustment curve giving the most probable heat-transfer law based on 
present information. 



Figure 11. - Theoretical streamlines for hori- 
zontal cylinder. The numbers "the values 

for '^/v • Scale is for Gr = 10^: For Gr = 10®, 
the distances from the surface are to be 
shortened by 1/3.16; for Gr = 10®, by l/lO. 


4 :^ 



Figure 12. - Theoretical isotherms for hori- 
zontal cylinder. The numbers are values of 
e/,0 . The scale is for Gr = 10^, to be 
shortened as in figure 11 for other values 
of Gr. 
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Figure 13. - Velocity profiles near cylinder referred to nondimens ional coordinates. 
Continuous curve represents theoretical solution. Points according to 
measurements of Jodlbauer for 2r = 9 centimeters; t^ = 99.2° C; t^^ = 18.1° C; 

Gr = 3.76x10^. 



Figure 14. - Temperature profiles near cylinder referred to nondimens ional 
coordinates. Continuous curve represents theoretical solution. Points 
according to measurements of Jodlbauer for 2r = 9 centimeters; ty = 99.2° C; 
t^= 18.1° C; Gr = 3.76x10^. 
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Figure 15. - Velocity profiles near cylinder. Continuous curve gives 

theoretical solution. Points according to measurements of Jodlbauer for 
2r = 5 centimeters; t^ = 104.6° C; too = 18.1° C; Gr = 6.54xlo5. 
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Figure 16. - Temperature profiles near cylinder. Continuous curve gives 
theoretical solution. Points according to measurements of Jodlbauer for 
2r = 5 centimeters; t^ = 104.6® C; too = 18.1° C; Gr = 6.54x10^. 
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Figure I 7 . - Azimuth, function g(x) giving variation of local heat-transfer 

coefficient or of a reciprocal distance from the wall which is characteristic 
for profile (for example^ boundary- layer thickness) along cylinder perimeter. 



Figure 18. - Comparison of heat transfer for rectangular plate and horizontal 
cylinder. For H = d^ under otherwise equal conditions^ the two diagrams 
give dimensional heat-transfer coefficients a to scale. 
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Figure 19. - Ratio of total heat transfer of horizontal cylinder and vertical 
plates (on both sides) of different heights under otherwise equal conditions. 
Ratio of boundary- layer Re at upper stagnation point of cylinder and upper 
edge of plate. Abscissa, heat transfer and Re; ordinate, plate height. 



Figure 20. - Vertical plate 100x100x1 centimeters for determining occurrence 
of turbulence for free convection. 
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•’igure 2i. - Horizontal cylinder (58.5-cm diam., 100-cm length) for determining 
occurrence of turbulence for free convection. 



Figure 22. - Schlieren photograph of heat transfer at the vertical plate. Time 
photograph (25 sec) with lens and camera. Surface temperature^ t^ = 100 C; 
tyy= 16.50 c* -b = y 48 millimeters; critical height (arrows), left 66.2 
centimeters, right 60 centimeters; critical Grashof number, 13.1x10®; 
critical boundary-layer Ee, 323. 
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Figure 23. - Schlieren photograph of the heat transfer at the horizontal cylinder. 
Time photographs (20 sec) with lens and camera, t^^ 102° C; too^ 18° C; 
h^ 752 millimeters; Gr of the cylinder^ 10.23:10°; critical azimuth (arrows) 
left 123.8°^ right 115.0°; critical boundary- layer Ee, 281. 
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Figure 24. - Schlieren photographs of heat transfer at one side of vertical 

plate. Simultaneous instantaneous photographs (about l/20 sec) with indirect 
illumination, t^^ 86° C; too 16.5^ C) b = 7^8 millimeters, (a) laminar; upper 
part of (b)^ transition; (c) and (d)^ turbulent. 
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Figure 25. - Schlieren photographs of heat transfer at the horizontal cylinder. 
Simultaneous instantaneous photographs (l/20 sec) with indirect illumination, 
ty = 102^ C; too= 18^ C; b = 752 millimeters; Gr^ 10.2x10®. (a) lower stag- 
nation pointy laminar; (b) azimuth of about 120^^ turbulent; (c) upper stagnation 
pointy upflow of warm air with current lead wires and so forth. 
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Figure 26. - Variation of position of start of turbulence with Gr 
number of cylinder, pomputed for critical boundary- layer Be, 285, 

/ 



Critical cylinder azimuth^ deg 


Figure 2’J . - Variation of position of start of turbulence Vith the cylinder 

diameter for surface temperature of 100*^ C in air at 20° C and 760 millimeters 
of mercury. Computed for critical boundary- layer Be^ 285. 
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